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ABSTRACT 


Described  is  an  experimental  technique  for  providing  long- 
duration  pulses  which  can  be  applied  to  a  part  or  all  of  a  simple 
structure  such  as  a  beam,  plate  or  cylinder.  The  technique  employs 
essentially  the  familiar  shock  tube  except  that  the  detonation  front  of 
a  gaseous  explosive  provides  the  shock  wave.  Many  pulse  shapes  can 
be  produced  by  placing  in  the  tube  and  against  the  target,  layers  of 
different  materials  such  as  styrofoam,  polyurethane,  and  Mylar, 
and  by  sending  the  explosively -induced  shock  wave  through  them  to 
the  target.  The  search  for  pulse  shapes  was  mainly  confined  to  those 
of  the  blast  type  but  the  generation  of  other  typeB  is  equally  feasible. 

Outlines  of  the  theoretical  treatments  of  four  problems  are 
given.  They  concern  the  responses  of  (a)  a  clamped  beam  to  a  blast 
pulse  uniformly  distributed  over  a  central  length,  (b)  a  simply  supported 
circular  plate  to  a  blast  pulse  uniformly  distributed  over  a  central 
circular  area,  (c)  a  clamped  circular  plate  to  a  rectangular  pulse  uni¬ 
formly  distributed  over  the  whole  plate,  and  (d)  a  clamped  circular 
plate  to  a  rectangular  pulse  uniformly  distributed  over  a  central  circular 
area.  Analytical  treatments  employ  the  rigid -plastic  theory  because  of 
interest  in  moderately  large  permanent  deformations  and  relative  sim¬ 
plicity  of  analysis. 

Long -duration  pulses  were  applied  to  clamped  beams  and  clamped 
circular  plates,  and  the  permanent  central  displacements  are  correlated 
with  theoretical  predictions. 

A  preliminary  «rudy  is  made  of  the  modification  of  the  applied 
pulse  due  to  the  mobility  of  the  target. 
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SECTION  I 


INTRODUCTION 

Re-entry  vehicles  are  liable  to  be  subjected  to  loads  ranging 
from  sharp  pulses,  which  for  structural  response  may  be  considered  as 
ideal  impulses,  to  pulses  with  durations  comparable  to  the  fundamental 
elastic  response  times.  Short-duration  loads  have  been  applied  in  tests 
on  actual  ICBM  structures  [  1.  1  -  1.4]  and  in  experimental  and  theo¬ 
retical  studies  of  simplified  ICBM  structural  components  such  as  cylin¬ 
ders,  plates,  and  beams  [  1.5  -  1.25].  As  a  result  much  is  known  about 
response  to  short  pulses.  Less  is  known  about  the  effects  of  long -duration 
pulses  [  1. 26  -  1. 29]  especially  when  these  act  only  on  parts  of  a  struc- 
ture.  One  reason  for  this  deficiency  is  the  dearth  of  meaningful  experi¬ 
mental  results  and  this  in  turn  iu  due  to  the  need  for  controlled  long- 
duration  loading  techniques. 

In  Section  II  is  described  such  a  technique  suitable  for  applying 
loads  to  part  of  the  surface  of  a  structure.  The  development  was 
guided  by  the  following  pulse  requirements: 

(1)  low  peak  pressures, 

(2)  durations  comparable  with  fundamental  elastic  response 
times, 

(3)  rise  time  much  smaller  than  duration, 

(4)  sharp  edges  to  the  pressure  distribution  whenever  appro¬ 
priate, 

(5)  pulse  shape  control,  and 

(6)  adaptability  for  use  over  larger  areas. 

The  long-duration  loading  technique  meets  these  requirements,  partly 
chosen  so  that  pulses  of  the  blast  type  (sudden  pressure  rise  and 
gradual  decay)  are  included,  but  in  addition  it  provides  medium  and  high 
peak  pressures,  medium -duration  pulses,  and  shapes  other  than  those 
associated  with  blast  pulses.  Section  II  also  includes  a  preliminary 


Numbers  in  brackets  indicate  references  listed  at  the  end  of  tnis  report. 
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study  of  the  effect  on  the  pulse  of  the  mobility  of  the  target  and  it  is 
concluded  that  the  effect  can  definitely  be  first  order. 

Section  III  is  a  theoretical  study  of  the  response  of  a  clamped 
beam  subjected  to  a  blast  pulse  uniformly  diet  .buted  over  a  central 
part  of  the  span.  On  the  basis  of  interest  in  moderately  large  perma¬ 
nent  deformations  the  rigid -plastic  theory  is  employed.  This  theory 
neglects  elastic  deformations  and  provides  a  relatively  simple  approach. 
The  clamped  beam  problem  was  chosen  primarily  because  of  its  sim¬ 
plicity  and  because  it  gives  much  insight  into  the  response  or  deforma¬ 
tion  mechanisms.  Some  structural  experiments  were  performed  on 
clamped  beams  but  since  the  reproducibility  required  some  improve¬ 
ment  the  results  are  regarded  as  preliminary.  They  indicate  that  the 
predicted  damage  is  from  2  to  4  times  the  actual  damage. 

Section  IV  is  a  theoretical  sutdy  of  the  response  of  a  simply 
supported  plate  subjected  to  a  blast  pulse  uniformly  distributed  over  a 
central  circular  area.  The  solution,  which  again  uses  the  rigid-plastic 
theory,  is  for  low  peak  pressures  only,  cr  rather  for  low  values  of  the 
ratio  of  the  peak  pressure  to  the  static  collapse  pressure  acting  on  the 
same  area. 

Section  V  is  a  theoretical  and  experimental  study  of  the  response 
of  a  clamped  circular  plate  subjected  to  a  blast  pulse  over  the  whole 
area  of  the  plate.  Using  rigid-plastic  theory  a  solution  is  presented  for 
the  special  case  of  a  rectangular  pulse.  In  the  experiments,  the  plates 
were  subjected  to  blast  pulses  having  an  exponential  decay.  Hence  true 
correlation  of  theoretical  and  experimental  damage  is  not  achieved.  It 
is  believed  that  the  rectangular  pulse  among  all  blast  pulses  with  the 
same  pressure  and  impulse  causes  the  greatest  damage  and  the  corre¬ 
lation  showed  that  the  predicted  damage  is  from  3  to  7  times  the  actual 
damage,  so  with  true  correlation  these  values  should  be  reduced.  This 
problem  was  chosen  as  a  preliminary  to  that  in  Section  VI,  where  the 
plate  is  partially  loaded,  because  of  the  complexity  of  these  kinds  of 
problems. 


Section  VI  is  a  theoretical  and  experimental  study  of  the  response 
01  a  clamped  circular  plate  subjected  to  a  blast  pulse  over  a  central 
circular  area  of  plate.  Using  rigid -plastic  theory  a  solution  is 
presented  for  the  special  case  of  a  rectangular  pulse.  In  the  experi¬ 
ments,  the  plates  were  subjected  to  pulses  which  may  be  considered 
approximately  rectangular  but  having  an  exponential  decay.  Correlation 
of  the  theoretical  and  experimental  central  deflections  shows  that  the 
predicted  damage  overestimates  over  the  wide  range  of  factors  from 
2  1/2  to  24.  The  large  factors  are  attributed  partly  to  the  fact  that 
the  experimental  pulse  is  not  rectangular  and  partly  to  the  use  of  thin 
plates  which,  due  to  their  mobility,  do  not  receive  the  full  impulse. 

Section  VII  is  a  description  of  a  simple  mathematical  model  of  the 
experimental  configuration  of  the  long-duration  loading  technique. 

The  model  successfully  describes  the  mechanics  of  the  operation  and 
forms  a  good  basis  for  further  refinement.  It  is  a  valuable  aid  to 
experimental  design  towards  achieving  the  pulses  desired. 

Section  VIII  is  a  study  of  the  effect  of  pulse  shape  on  simple  rigid - 
plastic  structures.  It  is  proved  that  for  a  certain  class  of  structures 
the  rectangular  pulse,  among  all  pulses  of  equal  peak  pressure  and 
impulse,  causes  the  greatest  damage.  It  is  also  indicated  that  the 
theorem  is  true  for  a  wider  class  of  rigid-plastic  structures.  The 
importance  of  this  study  lies  in  the  fact  that  it  is  much  easier  to 
analyze  structures  when  the  pulse  is  rectangular  and  the  result  serves 
as  an  upper  bound  on  the  damage  acquired. 
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SECTION  II 


EXPERIMENTS 


Long-Pulse  Technique 


In  developing  a  technique  for  providing  and  applying  long-duration 
pulses  the  following  properties  are  being  sought: 

(a)  pulse  duration  about  1/2  msec, 

(b)  short  rise  time  to  peak  pressure, 

(c)  monotonic  decay  following  peak  pressure, 

(d)  uniform  pressure  distribution  on  loaded  surface, 

(e)  each  load  unit  applicable  to  small  area  of  structure,  and 

(f)  units  capable  of  being  combined  for  larger  surface  loading. 


An  experimental  arrangement  which  provides  pulses  meeting  the 
above  requirements  when  they  are  applied  to  fixed  targets  is  shown  in 
Figs.  2.  1  and  2.2.  It  essentially  consists  of  a  shock  tube  in  which  the 
shock  is  the  detonation  front  of  a  50/50  gaseous  mixture  of  oxygen  and 
acetylene.  One  end  of  the  tube  is  placed  against  a  fixed  target  plate, 
the  junction  being  made  air-tight  by  means  of  a  rubber  gasket.  In  the 
fixed  endplate  is  mounted  a  pressure  transducer.  The  open  end  of  the 
tube  is  sealed  with  a  sheet  of  Mylar.  Inlet  and  outlet  hoses  for  the  gas 
mixture  pass  through  this  Mylar  sheet.  The  gas  is  detonated  at  the 
opposite  end  of  the  tube  to  the  pressure  transducer  [2.  l),  so  that  the 
detonation  front  travels  towards  the  endplate  producing  an  instantaneous 
rise  of  pressure  there  and  a  pulse  duration  depending  on  the  tube  length. 
Figure  2.  3  shows  a  typical  pulse  obtained  with  this  experimental  arrange¬ 
ment. 


One  limitation  of  this  arrangement  is  that  the  peak  pressure  is 
always  a  constant  for  a  given  gas  mixture  although  this  can  be  relaxed 
somewhat  by  varying  the  ratio  of  oxygen  to  acetylene.  Also  the  form 
of  the  decay  in  pressure  remains  the  same,  its  rapidity  depending  on 
the  tube  length.  Another  and  more  important  limitation  for  spot 
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FIG.  2.2  OXYACETYLENE  SHOCK  TUBE 


applications  is  that  separation  of  the  target  and  tube  allows  rarefaction 
waves  to  travel  over  the  target  surface  and  rapidly  curtail  the  pulse 
duration.  The  arrangement  is,  however,  suitable  for  clamped  plates 
where  the  plate  support  is  an  extension  of  the  shock  tube  (see  Section  V). 

An  arrangement  which  does  not  have  these  limitations  and  pro¬ 
vides  innumerable  pulse  shapes  consists  of  filling  the  end  of  the  tube 
next  to  the  target  with  various  materials  on  which  the  shock  wave 
from  the  detonated  gas  must  act.  One  such  arrangement  is  shown 
schematically  in  Fig.  2.4  where  alternate  layers  of  polyurethane  foam 
and  Mylar  form  a  spring-mass  system.  Between  the  gas  and  the  foam- 
Mylar  combination  is  a  layer  of  styrofoam  which  acts  as  a  light  piston 
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FIG.  2.3  ROD-GACC  RECORD  —  OXYACETYLENE  GAS  PULSE 


and  helps  to  spread  the  pressure  uniformly  over  the  face  of  the  target. 
By  changing  this  system  of  masses  and  nonlinear  springs  various  pulse 
shapes  can  be  obtained. 

It  is  desirable  that  the  compressed  spring -mass  system  should 
be  much  more  flexible  than  the  target  to  minimize  the  disturbance  of 
the  uniformity  of  the  pressure  distribution  as  the  target  deforms. 

This  property  influenced  the  choice  of  the  spring  material  and  mass 
geometry. 

In  Figs.  2.  5  to  2.  9  are  displayed  sequences  of  pulses  that  have 
been  obtained  by  varying  certain  parameters  of  the  spring-mass 
system  of  Fig.  2.4.  Attention  has  been  focussed  on  obtaining  pulses 
which  can  be  approximated  by  rectangles,  triangles,  exponential 
curves  or  simple  combinations  of  these.  Many  of  the  pressure  records 
display  the  oscillations  of  the  systems  but  they  are  of  such  high  fre¬ 
quency  that  a  structure  will  respond  only  to  the  mean  pressure. 


Details  of  the  pulses  and  config¬ 
urations  are  ihown  in  Tables  2.  1 
to  2.  5  on  the  pages  opposite  the 
pulse  records.  The  main  obser¬ 
vation  to  be  made  is  that  a  fairly 
high  degree  of  control  of  the  pulse 
shaping  can  be  obtained  with  the 
spring-mass  system.  Furthermore, 
the  pulse  shaping  experiments  are 
reproducible  and  provide  a  reasonably  uniform  pressure  distribution 
over  the  target  face.  This  last  point  was  established  by  using  three 
rod  gages,  located  at  the  center,  a  half-radius  point,  and  close  to  the 
edge  of  the  target  area. 

In  addition  to  the  cylindrical  tube  experiments,  spring-mass 
systems  were  tried  in  a  tube  of  rectangular  cross  section.  Similar 
pulse  shapes  were  obtained  out  both  the  reproducibility  and  uniformity 
of  pressure  distribution  were  found  to  be  unsatisfactory  (for  example, 
peak  pressure  variation  was  +  10%).  These  deficiencies  are  attributed 
to  the  cros 8 -sectional  dimensions  of  the  tube  used  in  the  experiments. 

A  rather  narrow  rectangle, 5  inches  by  1  inch,  was  used  resulting  in 
the  intrusion  of  the  edge  effects  (friction,  slight  lack  of  fit  of  spring- 
mass  system,  etc).  Should  a  rectangular  tube  with  the  same  aspect 
ratio  for  the  cross  section  be  required  for  an  experiment  it  is  suggested 
that  the  scale  be  increased  to  provide  10  inches  by  2  inches,  say.  Also, 
it  would  be  desirable  to  provide  finer  tolerances  by  using  machined 
parts.  Figure  2.  10  shows  the  experimental  arrangement  using  a  tube 
of  rectangular  cross  section  located  vertically  over  three  3-foot  rod 
gages. 
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FIG.  2.4  TYPICAL  SPRING-MASS  SYSTEM 


9 


Spring-Mass  System  for  Pulses  of  Figs.  2.  5(a)  to  (d) 
Configuration  Diagram 
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Table  2.  1  <—»  Experimental  Data 
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SCALES 
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Th*  Mylar  dleke  are  lO-mila  thick  and  ara  iocatad  at  l/2-inch  camera. 
Impulee  valuaa  ara  obtalnad  by  crude  curva  fitting. 
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lb/ftJ 

gm/cm 

Styrofoam 

4 

0.  064 

Polyurethane 

5 

o.  ogo 

Mylar 

<7.4 

1.4 

Observations 

1.  Trend  of  the  sequence  is  from  an  exponential  curve  to  a  triangle. 

2.  Pelk  pressures  increase  monotonicaiiy  from  about  1000  lb/in2 
to  1800  lb/in£, 

3.  Impulses  tend  to  decrease. 
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FIG.  2.5  ROD-GAGE  RECORDS 


Spring-Mass  System  for  Pulses  of  Figs.  ?..  6(a)  to  (f) 


Configuration  Diagram 


Table  2.2  — »  Experimental  Data 
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Th«  Mylar  dlaka  are  10-mlle  thick  id  art  locatad  at  l/2-lnch  cantara. 
'mpulai  valuaa  ara  obtained  by  crude  curva  fitting . 


Material 

Deaelty 

lb/ftJ 

gm/cm3 

Styrofoam 

4 

0.064 

Polyurethane 

5 

0.  080 

Mylar 

87.4 

1.4 

Observations 

1.  Tri.nd  of  the  sequence  is  from  a  slowly  decaying  exponential 
curve  toward  a  rectangle  followed  by  a  rapidly  decaying  exponen¬ 
tial  curve. 

2.  Peak  pressures  increase  monotonically  from  about  1000  lb/in2 
to  3200  lb /in2. 

3.  Impulses, gene  rally  higher  than  those  of  Fig.  2.  5, tend  to  increase. 
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FIG.  2.6  ROD-GAGE  RECORDS 
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Spring-Mass  System  for  Pulses  of  Figs.  2,  ?(a)  to  (d) 


Configuration  Diagram 


Table  2.  3  — —  Experimental  Data 
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Impaiaa  valuaa  ara  ebtalaad  by  credo  curva  fitting. 
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4 
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» 
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S7.  4 

1.4 

Observations 

1.  Trend  of  the  sequence  is  from  a  triangle  to  a  rectangle  followed 
by  a  triangle  or  exponential  curve. 
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2.  Peak  pressures  decrease  monotonically  from  about  1600  lb /in 
to  1200  lb/in2. 

3  Impulses  in  the  same  range  as  those  of  Fig.  2.  6  and  tend  to 
inc rease. 
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Observation! 
1.  Trend  < 


Trend  of  the  sequence  is  from  a  triangle  to  a  rectangle  followed 
by  an  exponential  curve. 


,2 

Peak  pressures  decrease  monotonically  from  about  750  lb/in 
to  500  lb/in2. 


Impulses  remain  fairly  steady  and  are  about  one -half  of  the  impulses 
in  Figs.  2.  5  to  2.  8  (only  about  half  the  length  of  gas  used). 
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FIG.  2.10  EXPERIMENTAL  ARRANGEMENT  —  SHOCK  TUBE 
OF  RECTANGULAR  CROSS  SECTION 

2.  Photodiode  Experiments 

In  order  to  assess  the  effect  on  the  pulses  of  the  mobility  of  the 
target,  pulses  from  the  fixed  target  configurations  associated  with 
Figs.  2.  8(d)  and  2.  9(e)  were  applied  to  free  disks  placed  against  the 
end  on  the  shock  tube  (Fig.  2.  11).  It  is  reasonable  to  suppose  that, 
at  least  in  the  early  stages  of  the  application  of  the  pressure,  the 
fixed  target  pressure  will  be  modified  most  if  the  plate  is  completely 
free  to  move  away. 
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CYLINDER 


LIGHT  SOURCE 


FIG.  2.11  DIAGRAM  OF  PHOTODIODE  SETUP  (during  motior) 

The  technique  consists  of  obtaining  an  x,  t  plot  for  the  disk  of 
sufficient  accuracy  to  allow  the  second  derivative  x  ,  or  the  accelera¬ 
tion,  to  be  computed  accurately.  Then,  by  Newton's  law  p(t)  =  mx 
it  is  possible  to  construct  the  pressure -tirr.e  relation.  Also,  for  a  pulse 
of  duration  T  ,  it  is  possible  to  find  the  impulse  since  I(T)  *  mx(T)  . 

A  diagram  of  the  experimental  setup  to  provide  an  accurate  x,  t 
plot  for  the  disk  is  shown  in  Fig.  2.  11.  A  lucite  "ladder"  in  a  light 
magnesium  frame  is  attached  to  the  center  of  the  disk  by  means  of  a 
ball  joint  (Fig.  2.  12)  and,  an  the  disk  r  ioves,  the  ladder  holder  is 
guided  along  teflon -lined  tracks  (Figs.  2.12,  2.  13  and  2.  14).  The 
strip  of  lucite  has  50  or  100  lines  per  inch  equally  spaced  on  one  surface, 
the  lines  being  perpendicular  to  the  direction  of  motion  and  having  a 
thickness  equal  to  the  space  between  them.  On  one  side  of  the  ladder 
is  located  a  powerful  light  source  while  on  the  opposite  side,  perpen¬ 
dicular  to  the  ladder,  is  a  lucite  rod.  The  curved  surface  of  the  rod 
is  coated  black  to  exclude  light  and  on  the  end  next  to  the  ladder  is  an 
array  of  lines  with  the  same  spacing  and  parallel  to  those  on  the  ladder. 
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FIG.  2.12  EXPERIMENTAL  ARRANGEMENT  FOR  PHOTODIODE  MEASUREMENTS 


PHOTO  OiOOE 


FIG.  2.14  EXPERIMENTAL  ARRANGEMENT  FOR  PHOTODIODE  MEASUREMENTS 


Each  line  is  completely  opaque  so  that  when  they  are  exactly  opposite 
each  other  the  maximum  amount  of  light  is  transmitted  along  the 
lucite  rod;  almost  complete  interference  or  darkness  prevails  when 
the  lines  are  exactly  opposite  the  spaces.  As  the  ladder  moves  across 
the  end  of  the  rod,  the  intensity  of  light  traveling  down  the  rod  varies 
and  is  converted  into  a  voltage  variation  by  means  of  a  photodiode 

located  at  the  other  end.  The  voltage  variation  is  recorded  on  oscillo¬ 
scopes. 

Figures  2.  15(a)  and  (b)  are  oscillograms  showing  the  voltage 
variation  when  a  ladder  having  50  lines  per  inch  is  attached  to  a 
1/2-inch-thick  aluminum  disk  which  receives  a  pulse  from  the  con¬ 
figuration  associated  with  the  record  ir  Fig.  2.  9(e).  The  sweep  rate 
is  20  nsec  per  Cm  and  by  using  appropriate  delays  in  two  dual  beam 
oscilloscopes  a  total  time  coverage  of  800  psec  is  achieved.  The 
distance  traveled  by  the  disk  during  the  time  interval  between  crests 
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FIG.  2.15  PHOTODIODE  RECORDS  EXPERIME  viTAL  CONFIGURATION  ASSOCIATED 
WITH  FIG.  2.9(c)  (Ladder  hat  50  linet/inch) 

or  troughs  is  0.  02  inch.  Thus  an  (x.t)  plot  for  the  disk  is  obtained 
and  is  shown  in  Fig.  2.  16  as  the  curve  labelled  "experiment." 
Additional  points  can  be  obtained  by  using  the  trace  between  the  crests 
and  troughs.  In  particular,  the  upper  trace  in  Fig.  2.  15(a)  gives  the 
initial  motion  of  the  disk. 

As  stated  earlier  the  second  derivative  of  the  (x,  t)  plot  gives 
the  pressure.  It  was  found  that  the  method  gave  the  initial  part  of 
the  pressure -time  diagram  satisfactorily,  but  around  the  peak  pressure 
and  at  later  times  it  appears  that  the  pressure  is  varying  quite  rapidly 
and  this  demands  greater  accuracy.  However,  the  ultimate  accuracy 
of  the  method  has  not  been  reached. 

As  an  inverse  method  one  can  take  the  pressure  pulse  acting 
on  the  fixed  target  (rod-gage  record),  idealize  it,  and  modify  it  so  that 
its  (x,t)  plot  falls  on  top  of  the  experimental  (x,  t)  curve.  In  Fig.  2.  16, 
curves  A  and  B  are  (x,  t)  plots  from  idealized  pulses  having  the  shapes 
shown  in  the  figure,  impulses  and  peak  pressures  equal  to  those  of  the 
fixed  target  pulse,  and  a  ramp  pressure  rise  taking  100  psec  to  reach 
peak  pressure.  Sue’  pulse  shapes  could  form  a  reasonable  starting 
point,  especially  for  early  times  where  the  mobility  of  the  target  has 
least  effect. 
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FIG.  2.16  (x,t)  PLOTS  FOR  FREE  DISK- 
CONFIGURATION  OF  FIG.  2.9(e) 


One  point  worthy  of  notice  is  that  the  impulse  imparted  to  the 
free  disk  is  considerably  less  than  that  on  a  fixed  target.  This  can  be 
seen  in  Fig.  2.  16  by  comparing  the  final  slope  of  the  curves  A  and  B 
with  the  final  average  slope  of  the  experimental  curve.  From  two 
photodiode  experiments  with  the  configuration  of  Fig.  2.  8(d)  the  final 
velocities  corresponded  to  impulses  of  0.208  and  0.219  lb-sec/in2  for 
a  fixed  target  impulse  of  0.  363  lb -sec /in2.  From  three  photodiode 
experiments  with  the  configuration  of  Fig.  2.9(e)  the  impulses  were 
0.  117  0.  119,  and  0.  120  lb-sec/in2  corresponding  to  the  fixed-target 
impulse  of  0.209  lb-sec/in2.  The  average  ratios  of  free-target  to 
fixed-target  impulses  are  respectively  0.  59  and  0.  57  which  represents 
a  considerable  reduction  of  impulse  due  to  full  mobility  of  the  target. 

3.  Mathematical  Model 

Each  assemblage  of  layered  media  used  to  sha  pulses  has  been 
regarded  as  and  is  called  a  spring-mass  system,  whereas  it  is  actu  y 
a  somewhat  more  complicated  system.  However,  a  mathematical  model 
waa  constructed  consisting  of  masses  and  linear  or  cubic  springs  which 
gives  a  reasonable  account  of  the  behavior  of  the  layered  media,  at  least 
for  the  few  cases  studied.  The  description  of  and  results  from  the 
mathematical  model  are  contained  in  Section  VII.  The  (x,  t)  plot  of  a 
freely  supported  disk  used  in  the  photodiode  experiment  with  the  con¬ 
figuration  of  Fig.  Z.  9(e)  is  shown  in  Fig.  2.  16. 
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4. 


Structural  Experiments 


Some  experinents  were  performed  in  each  of  which  a  central 
length  of  a  clamped  beam  was  subjected  to  loading  by  placing  the  beam 
against  the  target  end  of  the  shock  tube  of  rectangular  cross  section. 

The  pulse  applied  to  each  beam  was  taken  as  that  measured  against  a 
fixed  target  (rod-gage  record)  for  the  sam?  experimental  configuration. 

As  was  mentioned  above  at  the  end  of  Subsection  1,  the  pulses  obtained 
were  neither  sufficiently  reproducible  nor  uniform.  In  spite  of  this  a 
few  beams  of  6061 -T6  aluminum  were  loaded  to  demonstrate  the 
feasibility  of  the  structural  experiments.  The  results,  which  can  only 
be  regarded  as  preliminary,  indicate  that  the  ratio  of  the  theoretical 
to  experimental  final  central  deflections  range  from  about  3  to  4  for 
X  values  (ratio  of  peak  pressure  to  static  collapse  pressure)  ranging 
from  about  3  to  12,  the  theoretical  deflection  being  that  from  a  rec¬ 
tangular  pulse.  In  conclusion  it  should  be  stated  that  these  experiments 
could  almost  certainly  be  improved  by  using  wider  beams  and  hence  a 
wider  shock  tube.  This  remark  is  based  on  the  good  reproducibility 
and  uniformity  of  pressure  distribution  obtained  when  the  shock  tube  is 
cylindrical.  Experiment*  on  plates  are  described  at  the  ends  of  Sections 

V  and  VI. 

"j.  Large  Surface  Loading 

For  the  loading  of  larger  structural  surfaces  the  shock-tube  units, 
which  can  have  any  reasonable  cross  section,  can  be  placed  side  by 
side  to  cover  the  loading  area.  A  number  of  units  are  especially  re¬ 
quired  when  the  loading  surface  is  curved  to  ensure  that  the  wave 
strikes  the  surface  at  right  angles  and  to  ensure  that  the  layers  are 
not  buckled  or  crumpled  by  being  forced  to  occupy  a  smaller  area. 

Figure  2.  17  shows  a  section  through  a  suggested  two-unit  assembly  of 
shock  tubes  for  applying  a  load  to  a  cylindrical  surface.  Unlike  the  single 
units  the  assembly  shown  would  involve  a  detonation  wave  btriking  the 
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FIG.  2.17  SECTION  OF  A  TWO  UNIT 

ASSEMBLY  OF  SHOCK  TUBES 


styrofoam  piston  obliquely,  the  front  traveling  .'.long  the  top  of  the 
styrofoam  at  a  velocity  higher  than  the  detonation  velocity.  The 
pressure  pulse  for  this  setup  would  first  be  obtained  from  rod  gages 
and  possibly  photodiode  experiments  before  being  applied  to  the 
structure. 

It  should  be  possible  to  approximate  a  loading  which  varies  not 
only  with  time  but  with  position  on  the  loaded  surface  by  using  an 
assemblage  of  shock  tubes  each  with  its  own  spring-mass  configuration. 
A  cosine  distribution  of  loading  around  one -half  of  a  cylinder  appears 
quite  feasible. 
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SECTION  III 


CLAMPED  BEAM 


1 .  Introduction 

The  problem  treated  here  is  the  response  of  a  clamped  beam  of 
rigid-plastic  material  subjected  along  a  central  portion  of  its  span  to 
a  pressure  which  is  constant  along  the  beam  but  varying  with  time  in 
the  form  of  a  blast  pulse.  Figure  3.  1  illustrates  the  problem.  A 
blast  pulse  is  here  defined  as  a  pulse  with  a  time-dependent,  pressure 
satisfying 

t 

tp(t)  S  J  p(T)dT  (3.  1) 

o 


FIG.  3.1  CLAMPED  BEAM  PROBLEM 


The  main  characteristics  of  such  a  pulse  are  the  instantaneous 
rise  to  peak  pressure  and  a  decay  in  accordance  with  condition  (3.  1). 

When  the  moment-curvature  relation  for  a  beam  is  approximated 
by  that  corresponding  to  a  beam  made  of  a  rigid-ideally  plastic  material, 
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no  deformations  and  hence  no  curvature  is  possible  until  the  bending 
moment  at  a  cross  section  is  equal  to  the  fully  plastic  moment,  where¬ 
upon  the  curvature  becomes  unbounded.  This  property  determines 
failure  in  rigid -plastic  structures  by  means  of  mechanisms. 

In  the  Semiannual  Technical  Report  No.  1  it  was  shown  that  the 
clamped  beam  can  fail  in  one  of  four  types  of  mechanisms.  The  param¬ 
eters  determining  the  initial  mechanism  are  the  peak  pressure  of  the 
pulse  and  the  loaded  length  of  the  beam.  In  the  following  section  the 
equations  of  motion  ant’  their  solutions  are  presented  for  each  mecha¬ 
nism.  Two  of  these  mechanisms  are  the  subject  of  a  oaper  by  M. 
Conroy  [3.  1 } .  Although  the  work  primarily  concerns  oeams  of  infinite 
length,  its  application  for  a  clamped  beam  can  easily  be  made. 

Section  III  considers  the  deformations  caused  by  three  types 
of  blast  pulses:  rectangular,  triangular  and  exponential.  A  derivation 
of  formulas  for  the  deformed  shape  of  the  beam  is  presented  only  for 
the  case  of  a  rectangular  pulse.  Formulas  for  the  midspan  deflections 
are  contained  in  Appendix  A.  For  the  rectangular  pulse  only,  formulas 
for  the  deformed  shape  are  contained  in  Appendix  B.  The  variation  of 
the  midspan  deflection  with  peak  pressure,  impulse  and  loaded  length 
for  all  three  pulses  are  shown  in  Figs.  3.  3  to  3.  1 1. 

2.  Mechanisms  of  Deformation 

(a)  Mechanism  1 


Peak  pressure!!,  pm  ,  which  are  slightly  above  the  static 
collapse  pressure  p(  ,  cause  small  inertia  forces  and  the  beam  deforms 
in  the  static  collapse  mechanism.  The  cross  sections  of  the  beam  at 
the  supports  and  at  midspan  carry  a  fully  plastic  moment,  M(>  , 
and  the  resulting  motion  of  each  half-beam  is  a  rotation  as  a  rigid  body 
about  the  supports.  The  corresponding  velocity  distribution  along  a 
half-beam  ,  0  s  x  s  l  ,  is 

y  (x,  t)  *  ju(tMt"x)  (3.2) 
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ft* 


*  ✓ 

lip 


where  U)(t)  is  the  angular  velocity  about  the  ends, 
motion  can  readily  be  written  in  the  form 


The  equation  of 


uu  = 


6M  p 
os 


mi 


?  P 


(  3-  -  1  ) 

s 


(3.3) 


where  p#  is  the  static  collapse  pressure  given  by 

4M 


*  S0«2-S0> 


(3.4) 


and  m  is  the  mass  per  unit  length  of  beam. 


The  equilibrium  equation  for  a  length  x  of  the  half-beam 
(0  s  x  £  i  )  i8 

X  X. 


1 

M  =  Mo  ~J  J  (p-my)dx'  dx 


(3.5) 


When  y  .  found  by  differentiating  (3.2)  and  using  the 
result  (3.3)  for  i  ,  is  substituted  in  (3.5)  and  the  integration  performed, 
(3.5)  provides  the  following  relations  for  the  tending  moment 


=  1  +  rf-PS2  -  §3)  (-^- 


■  -  S 


M  .  , 

K  '  1  p 


+  ^(3?2  -53)  (^i-) 


0  5  5  S’  §, 


50  s  ^  s  1 


(3.6) 


The  conditions  which  ensure  no  violation  of  the  yield 


condition,  |  M|  £  Mq  ,  are  now  stated. 


From  expressions  (3,  6)  it  can  be  shown  that  a  necessary 
and  sufficient  condition  for  a  maximum  moment  at  midspan  is 


s 


3*o  <2  -  ?o> 


'O 


p  3S0<2-50> 


s  i  - 1 /yr=  .433 

(3.7) 


§0  *  i  -  i/yr 
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For  the  case  wuere  the  loading  acts  over  more  than  half 
the  span  and  X  >  Xj  ■  the  yield  condition  is  violated  in  the  neighborhood 
of  5  =  0.  This  suggests  a  mechanism  of  deformation  which  retains 
the  plastic  hinges  at  the  supports  but  has  two  plastic  hinges  which 
travel  towards  the  center.  The  decreasing  central  section  of  the  beam 
will  translate  downward  while  the  outer  portions  rotate  as  rigid  bodies 
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about  the  supports.  After  the  moving  hinges  meet  at  the  center  the 
beam  will  deform  according  to  mechanism  1. 


Denoting  the  position  of  the  moving  hinge  by  x  (t)  ,  the 
velocity  distribution  meeting  the  above  description  is 


y(x,  t) 


yG(t) 


y  (t)  (-£-li) 

1  o'  'x  -l' 
o 


0  £  X  £ 


X  £  X  £  l 

o 


(3.  10) 


The  equations  of  motion  for  the  half  beam  are 


mV0  =  P  (3.11) 

-  xq)2  tv  =  £(a  -  xq)  [Zl  -  a  -  xq)  -  2Mq  (3.  12) 

where  id  =  YQl(l  -  xq)  is  the  angular  velocity  of  the  outer  portions 
of  the  beam. 


where 


By  substituting  the  velocity  field  (3.  10)  in  the  moment 
expression  (3.  5)  and  carrying  cut  the  integration  the  moment  distri¬ 
bution  is  found  to  be 
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•  4  p"  f** 


Mo  '  £  (x‘xo)2  -  * 

Mo‘f  (a'Xo)(2x‘a‘xo)’IFcP 


0  s  X  s  X 


x  s;  x  u  a 
o 


a  sx  s  ( 


(3.  15) 


where 


3  2  3  2 

Co=x  -3x  (  -  3xx  +  6xx  (,+  2x  -  3x  l 

o  o  o  o 


Using  (3.  15)  it  can  be  shown  that  the  yield  condition  will 
not  be  violated  so  long  as  the  pulse  obeys  the  condition  (3.  1)  and 
X  =  p  /p  is  bounded  by  the  values  X.  and  X,  ;  that  is 

Tt\  8  1  b 


where 


Xj  *  X  £  X2 


3?o<2-  V 
?  -  - 2 — 


S  *  1/2 

o 


and  Xj  is  given  by  (3.  9) . 


(c)  Mechanism  2b  ?  s  1/2 

o 


When  the  loaded  length  is  less  than  the  half  span  and 
X  >  Xj  the  yield  condition  is  violated  at  the  supports.  A  mechanism 
of  deformation  is  thus  considered  which  jias  a  plastic  hinge  at  mid¬ 
span  and  two  hinges  which  travel  towards  the  supports.  The  only 
motion  is  that  of  the  inner  portion  of  the  half  beam,  0  s  x  £  x^  , 
which  rotates  as  a  rigid  body  about  the  moving  hinge  at  x  =  x^(t)  . 
After  the  moving  hinges  reach  the  supports  the  structure  will  deform 
according  to  mechanism  1. 


34 
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In  accordance  with  the  above  motion,  the  velocity  distri¬ 


bution  is 


y  (x,t)  = 


y0(‘)(l-^-)  Osxsx, 


Xj  S  X  s  l 


where  yQ  is  the  velocity  at  the  center  of  the  beam. 

The  momentum  equations  for  the  half -beam  are 


(3.  17 


m  dt(xlyo>  =  2ap 


(3.  18 


d  .  2  •  ,  ,2 

m  "ar  (xi  y0)  =  3a  p  +  12MC 


(3.  19 


With  the  initial  conditions  of  zero  displacement  and  velocity 
the  solution  of  (3.  18)  and  (3.  19)  is 


4a2  I2 (t) 

-  2 - 

12MQmt  +  3a  ml(t) 


(3.20 


x  .  3.  6Mo‘ 

xi  ■  ~r  +  nr 


(3. 21; 


This  mechanism  was  investigated  in  [3.  l]  and  it  was  found 
that  the  yield  condition  will  not  be  violated  ro  long  as  (3.  1)  is  satisfied 


a 


(3.22] 


Hence 


In  (3.22),  I (t ) / 1 — ►  p  as  t — ►  0,  so  that  p  s  1»’M  / a 

m  1  m  o 


where 


X2  =  3(  ?, 


Sn  *  1/2 

o 


(3.23) 


and  X  j  io  given  by  (3.  9) . 

(d)  Mechanism  3 

For  sufficiently  high  pressures  such  that  >  >  X^,  * 
mechanism  2  violates  the  yield  condition  either  at  the  supports  (if 
f  £  1/2)  or  at  mid-span  (if  ?  2  1/2).  This  suggests  a  mechanism 

of  deformation  in  which  moving  plastic  hinges  are  formed  on  both 
sides  of  the  loading  boundary,  x  =  a  .  The  hinge  at  x  =  xq  and  under 
the  loaded  segment  of  the  beam  moves  towards  midspan  while  the 
other  at  x  =  x^  moves  towards  the  support.  If  the  support  receives 
the  approaching  hinge  first,  the  structure  enters  mechanism  2a, 
whereas  if  the  mid-span  receives  the  other  approaching  hinge  first, 
the  structure  enters  mechanism  2b. 


The  motion  of  the  half  beam  is  as  follows:  the  inner 

segment,  0  £  x  £  xq  translates  downward  while  the  midsection, 

x  £  x  £  x,  ,  rotates  about  the  outer  hinge,  x.  .  The  velocity  field 
o  1  1 

can  therefore  be  written  as 

0  £  x  £  x 

o 

x  £  x  £  x,  (3. 24) 

O  1 

/ 

Xj  *  x  £ l 


y(x,  t)  =  { 


X  .  -  X 

y0  (r~^r) 

°  51 1  o 


The  momentum  equations  for  the  half- beam  are 
t  l 


JaplTMT  =  j  my  dx 


and 


k 

J  (2MQ  +  p^-)dT  =  J  mxydx 
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with  the  solution 


m  y 
'  o 

=  J  I(T)d7 

(3.25) 

0 

X  = 

0 

12M  t  ,  ,, 
a  -  (—jJL)1/* 

(3.26) 

X1  = 

12M  t  ,  ,, 

•  ti-fi1'2 

(3.27) 

In  [  3.  l]  it  is  shown  that  the  yield  condition  will  not  be 
violated  if  the  loading  is  a  blast  pulse  satisfying  (3.  1)  and 


X  ^  3(  — ) 

*  A 

1/2 

u 

(3.28) 

n  3'o<2-?o> 

S  *  1/2 

(1-?/ 

o 

To  find  the  initial  locations  of  x^ 

and  x_  set  I  =  p  t 

2  m 

26)  and  (3.27)  which  is  approximately  true  for  small  t  .  Then 

•V°>  -  *0  + 

*1<°>  =  *o  •  <UMo'Pm>‘/2 


Hence  as  \m  or  pm  increases  indefinitely  the  initial 
positions  of  Xj  and  move  closer  to  xq  ,  coinciding  with  xq 
for  an  ideal  impulse. 

(e)  Summary  of  Mechanisms 

The  results  obtained  for  the  modes  of  failure  giving 
dependence  upon  dimensionless  peak  pressure  and  loaded  length  are 


shown  in  Fig.  3.2,  For  convenience,  and  are  given  below. 


350(2  -?0)/[3?0{2-?0)-2j 

3(2  -  ?Q)/t  3(2  -  5C)  -  4] 


3?0(2-;o,/l-50)2 

3<2-5o'?o 


1/2  s  5  si 

0  <  ?  s  1/2 
’o 

1/2  s  f  si 
o 

0  <  §  s  1/2 


Once  a  blast  pulse  is  known  its  peak  pressure  and  the  loaded 
length  of  beam  can  be  plotted  on  Fig.  3.2  to  give  the  initial  mechanism. 
Knowing  this,  th?  analysis  can  proceed  directly  towards  obtaining 
permanent  deformations. 

3.  Permanent  Deformations 

(a)  Introduction 

In  this  Section  the  permanent  deformations  are  considered 
for  three  blast  pulses:  rectangular,  triangular  and  exponential.  These 
particular  pulses  need  only  two  parameters  for  their  complete 
description.  The  parameters  chosen  are  the  peak  pressure,  pm  , 
and  the  total  impulse  per  unit  beam  length,  Ij  ,  defined  as 

? 

1\  -■  J  P(t)dt 

o 

Since  the  method  of  solution  is  identical  for  all  three 
loadings,  expressions  are  derived  for  the  final  deformed  shape  due  to 
a  rectangular  pulse  causing  initial  deformation  by  mechanism  2a.  For 
the  other  initial  mechanisms  only  the  results  are  given  and  they  may  be 
found  in  Appendix  B.  The  central  deflection  formulas  for  all  three 
pulses  are  contained  in  Appendix  A. 
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(b)  Rectangular  Pulse 


The  loading  under  consideration  is  given  by 

Cpm  0  s  t  s  t 

J  m  o 

P=lo 

where 

t  =  I  /p 

o  1  m 

A  special  feature  associated  with  a  rectangular  pulse  is 
that  all  plastic  hinges  are  stationary  while  it  is  acting  on  the  beam. 

This  simplifies  the  analysis  since  the  quantity  I(t)  ,  appearing  in  the 
equations  of  Section  III-2,  canbe  replaced  by  the  constant  1^  , 

The  derivation  of  expressions  for  the  final  deformed  shape 
is  confined,  for  simplicity,  to  the  initial  mechanism  2a.  For  the 
remaining  initial  mechanisms  the  expressions  are  listed  in  Appendix  B. 

Mechanism  2a  :  X  .  <  X  s  X  1/2  t  f  *  1 

The  expressions  for  the  final  shape  of  the  half  beam  will 
be  quite  different  on  either  side  of  the  initial  hinge  position  x  =  x  (0) 
for  two  reasons,.  Firstly,  a  discontinuity  of  slope  at  x  =  x  (0)  is 
created  during  phase  1  because  a  rectangular  pulse  causes  the  hinge 
to  remain  stationary.  Secondly,  during  phase  2  the  hinge  initially  at 
x  =  xo(0)  travels  through  the  region  0  s  x  s  xq(0)  toward  the  center 
of  the  beam. 

During  phase  1,  the  equation  of  motion  in  the  region 
0  s  x  s  xo(0)  is  my  =  p^  which  upon  integrating  twice  between  t  =  0 

and  t  =  t  =  I/p  yields 
o  m 

y(x,  tQ)  =  I“/2mpm  0  s  x  s  xq(0)  (3.29) 
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The  portion  of  beam  xq{0)  *  x  s  l  rotates  about  the  support  aa  a 
rigid  body.  Hence 

y(x*tQ)  =  (I^/2mpm){{,-x)/U-x  (0)]  x  (0)  sxsi  (3.30) 


During  phase  2  let  the  time  be  t  when  the  moving  hinge 
reaches  t  point  in  the  segment  0  i  x  s  xq(0)  .  The  point  moves  at 
a  constant  velocity  I/m  during  the  time  t  -t  after  which  it  moves 
as  part  of  the  rigid  portion  rotating  about  the  support  and  ior  a  time 
*1  T  *  *  1  denoting  the  end  of  phase  2.  Now  the  equations  of  motion 
for  the  parts  on  either  side  of  x  *  xq  are  yQ  =  l/m  and  m({,-x  )3a./3  = 
-2Mq  from  which  xq  *  -6Mq/I«, -xq)  and  ft-x)2  -  U-xo(0)]2° 

6Mq(t  -  tQ)/I  .  The  velocity  of  points  on  the  rigid  rotating  portion  is 
y=  I(f,-x)/m({,-xo).  Hence  the  deflections  which  occur  during  phase  2 


are,  by  integration  of  velocities, 

I, 

y(x.t.)  -y(x.t  )  -  -Mt-tlt 


i.  i 

— [T-to)  +  (t-x)  J  - 

J  > 


»  X  (t-X  ) 

x  o’  O  ' 


=  I^[R-x)2-  t-x  (0)  +x(t-x)]/6mM 


y(x,tj)  -  y(x, t  ) 


0  i  x  sxQ(0)  (3.31) 


T^T  =  ^  "X)  f  • 

"  o  r  x 

xJO) 


xo(t-xo> 


llxo<0)  •  (£-x)/6mM  x  (0)  s  x  s,  i  (3.32) 
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Figures  3.4(a)  to  (d)  show  the  same  curves  as  those  of  Figs. 

3.  3(a)  to  (c)  only  they  have  been  regrouped.  For  a  given  loaded  length 
the  X  versus  v  curve  for  each  of  the  three  types  of  pulse  are  shown 
side  by  side.  Again  regarding  the  impulse  as  the  same  for  each  type 
of  pulse,  Figs.  3.  4(a)  to  (d)  show  that  for  a  given  value  of  X.  the 
central  deflections  are  greatest  when  the  pulse  shape  rectangular 
and  least  when  it  is  exponential.  For  low  values  of  X  the  spread  of 
6  is  quite  large  whereas,  as  expected,  for  high  values  of  X  the  spread 
is  small.  The  latter  observation  is  of  course  due  to  each  pulse  tending 
toward  an  ideal  impulse. 

Figures  3.  5(a)  to  (c)  give  the  same  information  as  that  in 
Figs.  3.3(a)  to  (c)  but  in  a  form  convenient  for  studying  the  variation 
of  X  with  the  impulse  .  They  are  essentially  pressure -impulse 
diagrams  (p-I  diagrams), each  curve  showing  how  the  pressure  and 
impulse  must  vary  to  achieve  a  given  central  deflection  6  .  It  can 
be  seen  that  above  certain  values  of  X  the  peak  pressure  can  be  varied 
significantly  with  very  little  change  of  impulse  required  to  maintain  6 
(see  Conclusion  4). 


Figures  3.  6(a)  to  (c)  show  the  same  curves  as  those  of  Figs.  3.  5(a) 
to  (b)  only  they  are  regrouped  according  to  loaded  length  so  that  the 
effect  of  the  pulse  shape  may  be  seen  more  clearly.  For  any  fixed 
peak  pressure,  or  X  value,  the  impulse  required  to  produce  \  given 
6  is  least  when  the  pulse  is  rectangular  and  greatest  when  it  is  expo¬ 
nential.  The  differences  are  most  pronounced  when  X  is  low,  especially 
in  the  range  1  <  X  <  2  . 

Figures  3.  7(a)  to  (c)  are  essentially  another  form  of  pressure- 
impulse  diagram  (p-I  diagram).  The  impulse  Ij  has  been  rendered 
dimensionless  by  dividing  by  I.  ,  which  is  the  ideal  impulse  required 
to  give  the  same  central  deflection  as  .  Consequently  all  the  curves 
have  a  vertical  asymptote  through  Ij/I.  =  1.  For  any  given  central 
deflection  6  the  curves  show  the  relationship  between  X  and  Ij/I.  . 

The  curves  are  of  course,  similar  to  those  of  Figs.  3.  5(a)  to  (c). 
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Figures  3.8(a)  to  (c)  show  the  same  curves  as  those  in  Figs. 

3.  7(a)  to  (c)  but  they  are  regrouped  according  to  loaded  length  to 
bring  out  the  effect  of  the  pulse  shapes.  They  arc  similar  to  those 
of  Figs.  3.  6(a)  to  (c). 

Figure  3.  9  shows  the  relationship  between  the  peak  pressure 

and  the  central  deflection  by  means  of  the  parameters  p  and  v  . 

This  is  an  alternate  way  of  representing  the  Information  in  Figs.  3.  3(a) 

to  (c)  to  give  the  peak  pressure  directly  instead  of  from  X.  =  p  /p  . 

m  s 

Figures  3.  10(a)  to  (c)  show  the  relationship  between  the  central 
deflection  6  and  the  loaded  length  a  by  means  of  the  parameters  v 
and  5q  with  each  curve  representing  a  constant  peak  pressure  charac¬ 
terized  by  the  parameter  p  .  A  vertical  line  drawn  through  some 
chosen  value  of  5^  intersects  the  p  curve  for  the  peak  pressure  of 
interest.  Then,  for  the  impulse  of  interest  the  central  deflection  can 
be  calculated  from  the  ordinate  value  of  v  .  Figures  3.  10  can  there¬ 
fore  be  regarded  as  design  curves  for  clamped  beams  of  rigid-plastic 
material.  Where  the  p  curves  intersect  the  horizontal  axis  v  =  0 
gives  an  idea  of  the  loaded  length  below  which  no  deflection  occurs. 

For  example,  if  p  =  10  more  than  one -fifth  of  the  span  must  be  loaded 
to  give  a  permanent  central  deflection.  Now  in  this  region  elastic  effects 
become  important  and  so  it  is  probable  that  one -fifth  is  a  lower  bound 
of  the  fraction  of  span  that  must  be  loaded  to  give  a  permanent  deflection  . 
If  a  beam  is  subjected  to  a  pulse  with  a  shape  that  can  be  approximated 
by  a  rectangular,  triangular  or  exponential  shape  of  the  same  impulse 
and  peak  pressure  then,  for  each  impulse  and  central  deflection  ,  v  is 
determined.  A  horizontal  line  through  this  value  of  v  cuts  p  curves 
which  give  the  relationship  between  the  peak  pressure  and  loaded  length 
to  maintain  the  central  deflection  (regard  p  curves  as  contours).  The 
loaded  length  must  always  be  greater  than  the  value  at  the  intersection 
of  the  horizontal  line  through  v  with  the  curve  p  =  ©  which  represents 
the  ideal  impulse  case. 


Figure  3.  11  shows  the  curves  p  =  20  for  each  type  of  pulse 
plotted  in  the  v,  plane  and  shows,  for  a  given  ?  ,  how  the  shape 

effects  the  final  deflection  when  the  peak  pressure  and  impulse  remain 
fixed. 


5.  Conclusions 


From  the  information  given  in  Figs.  3.  3  to  3.  11  the  following 
conclusions  can  be  drawn. 

(1)  For  a  given  impulse  Ij  and  loaded  length  50  ,  the  central 
deflection  increase  j  monotonically  with  pressure  pm  (or  \), 
and  becomes  t.  maximum  when  the  pressure  is  infinite,  that 
is,  when  the  given  impulse  is  applied  as  an  ideal  impulse. 
This  can  be  seen  best  in  Figs.  3.  *5  to  3.  8. 

(2)  For  a  given  impulse  Ij  ,  Table  3.  1  shows  approximate 

values  of  \  ,  corresponding  to  pulse  shape  and  loaded 
length  ?0  ,  above  which  over  80  percent  of  the  maximum 
central  deflection  (the  deflection  6 ,  when  Ij  is  ideal, 
i.e.,  1 1  =  Ii )  is  obtained.  Below  tne  listed  values  of  \ 

the  decrease  of  6  with  \  is  quite  pronounced.  This 
behavior  can  be  observed  best  in  Figa.  3.  3  and  3.  4.  In 
Table  3.  1  scmt  lower  bound  values  of  \i  are  also  listed. 


Table  3.  1 

LOWER  BOUNDS  FOR  \  AND  li  GIVING  6  >  0.  8  6. 
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(3)  The  effec  of  pulse  shapes  on  the  midspan  deflection  6 
when  the  peak  pressure  pm  ,  impulse  Ij  and  the  loaded 
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length  ratio  5  are  given  can  be  seen  in  Figs.  3.  4  and  3.  9. 
Rectangular  pvftses  cause  the  greatest  midspan  deflection 
and  the  exponential  pulsej  cause  the  least. 

(4)  For  a  given  midspan  deflection  6,  as  the  pressure  ratio 
X  decreases  from  infinity  of  some  value  X jj,  the  deflection 
6  is  maintained  with  less  than  a  10  percent  increase  in 
impulse  over  the  ideal  impulse  (1  <  Ij  <  1.1  for  X jj  <  X  <  00). 
Values  of  XL  obtainable  from  Fig.  S.  7  or  3.8  are  shown 
in  Table  3.2  for  the  three  types  of  pulses  and  for  %Q-  1/4, 
1/2,  and  1  .  For  values  of  X  below  XL  a  significant 
increase  in  impulse  is  required  to  maintain  6,  especially 
in  the  range  1  <  X  <  2  . 

Table  3.2 


LOWER  BOUNDS  OF  X  REQUIRING 
Ij/Ij  <  1.1  TO  MAINTAIN  6 
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FIG.  3.3(o)  PRESSURE  RATIO  vs.  MIDSPAN  DEFORMATIONS:  RECTANGULAR  PULSE 
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FIG.  3.3(c)  PRESSURE  RATIO  ...  MIDSPAN  DEFORMATIONS:  EXPONENTIAL  PULSE 


9 


6«-«*4W 


FIG.  3.4(a)  COMPARISON  OF  PRESSURE  RATIO  v«.  MIDSPAN  DEFORMATIONS 
FOR  RECTANGULAR,  TRIANGULAR,  AND  EXPONENTIAL  PULSES 
Loading  Ungth  ratio:  one  quarter 


FIG.  3.4(b)  COMPARISON  OF  PRESSURE  RATIO  v».  MIDSPAN  DEFORMATIONS 
FOR  RECTANGULAR,  TRIANGULAR,  AND  EXPONENTIAL  PULSES 
Loading  Ungth  ratio:  on#  half 


FIG.  3.4(c)  COMPARISON  FOR  PRESSURE  RATIO  v*.  MIDSPAN  DEFORMATIONS 
FOR  RECTANGULAR,  TRIANGULAR.  AND  EXPONENTIAL  PULSES 

Loading  length  ratio:  three  quartan 


0  004  QM  ait  OH  020 


FIG.  3.4(d)  COMPARISON  FOR  PRESSURE  RATIO  v*.  MIDSPAN  DEFORMATIONS 
FOR  RECTANGULAR,  TRIANGULAR,  AND  EXPONENTIAL  PULSES 
Loading  length  ratio:  total  span 
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FIG.  3.5(b)  PRESSURE  RATIO  v».  IMPULSE:  TRIANGULAR  PULSE 
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FIG.  3.5(c)  PRESSURE  RATIO  v».  IMPULSE:  EXPONENTIAL  PULSE 


FIG.  3.6(o)  PRESSURE  RATIO  v«.  IMPULSE:  ALL  PULSES 


FIG.  3.6(b)  PRESSURE  RATIO  vs.  IMPULSE:  ALL  PULSES  £0  -  K 
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FIG.  3.7( 
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)IA3RAM:  RECTANGULAR  PULSE 

FIG.  3.7(b)  PRESSURE-IMPULSE  DIAGRAM:  TRIANGULAR  PULSE 
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FIG.  3.7(c)  PRESSURE-IMPULSE  DIAGRAM:  EXPONENTIAL  PULSE 


3.8(a)  PRESSURE-IMPULSE  DIAGRAM:  ALL  PULSES  <£, 
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FIG.  3.8(b) 


PRESSURE-IMPULSE  DIAGRAM:  ALL  PULSES  £ 0 
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FIG.  3.8(c)  PRESSURE-IMPULSE  DIAGRAM:  ALL  PULSES  £, 


FIG.  3.10(b)  DEFLECTION,  LOADED 
TRIANGULAR  PULSE 


FIG.  3.10(c)  DEFLECTION,  LOADED 
EXPONENTIAL  PULSE 


FIG.  3.11  DEFLECTION,  LOADED  LENGTH,  PRESSURE  RELATIONSHIP: 
ALL  PULSES  fi  m  20 


NOMENCLATURE 


load  half  length 
impulse  at  time  t  <  t 
total  impulse 
ideal  impulse 
half  span 

mass  per  unit  length 
bending  moment 
fully  plastic  moment 


SECTION  IV 


SIMPLY  SUPPORTED  RIGID-PLASTIC  PLATES 
UNDER  BLAST  LOADING 


1 .  Introduction 

A  simply  supported  circular  plate  of  rigid-plastic  material 
subjected  to  a  blast  load  of  low  peak  pressure  acting  on  a  central 
circular  ar-a  is  analyzed  in  thi*  section  (Fig-  4.  1).  The  pressure 
is  assumed  to  rise  instantaneously  to  its  peak  pressure  p^  and  there¬ 
after  to  decay  monotonically.  Only  one  deformation  mechanism  is 
investigated,  namely,  that  corresponding  to  static  collapse.  Conse¬ 
quently,  only  the  response  of  the  plate  to  low  peak  pressures  is 
presented  below. 

Yielding  is  assumed  to  occur  in  bending  according  to  the  Tresca 
yield  condition  and  the  associated  flow  rule  (Fig.  4.2)  [4.  l],  membrane 
action  being  neglected.  In  Fig.  4.  2,  and  Mg  are  the  radial  and 

circumferential  components  of  bending  moment.  ky  and  kg  are  the 
corresponding  components  of  rate  of  curvature.  Positive  bending 
moments  and  shear  force  Q  are  shown  in  Fig.  4.  3. 


FIG.  4.1  CIRCULAR  PLATE 
PROBLEM 


FIG.  4.2  TRESCA  YIELD 
HEXAGON 


FIG.  4.3  PLATE  ELEMENT 


2.  Deformation  Under  Low  Pea?T  Prer  lures 

The  deformation  mechanism  is  assumed  that  of  static  collapse, 
namely, 

w(r,t)  =  W(t)(l-r/R)  (4.1) 

where  R  is  the  plate  radius.  Deflection  formula  (4.  1)  describes  a 
deformation  in  which,  at  any  instant,  the  plate  is  a  shallow  cone. 
According  to  (4.  1)  the  rates  of  curvature  are 

H  =  -  w  *  0  and  kd  =  -  w  /r  =  W/Rr  (4.2) 

r  rr  0  r 

The  flow  rule  as  described  by  Fig.  4.  2  and  results  (4.  2)  dictate 

that  the  plate  is  in  the  plastic  regime  AB.  In  fact,  at  the  plate  center 

M0  =  M  =  M  which  corresponds  to  the  regime  A;  at  the  support  M  =0 
u  r  o  t 

which  corresponds  to  regime  B.  Thus  the  bending  moments  corres¬ 
ponding  to  Eq.  (4.2)  are 

Mg  =  Mo  and  0  *  s  MQ  (4.3) 

where  MQ  i»  the  fully  plastic  moment  per  unit  length  of  cross  section. 

With  the  aid  of  Fig.  4.  1  the  equation  of  motion  can  be  derived  in 
the  form 

r 

M0-^(rMr)  -y  (p-mw)rdr  (4.4) 

o 

where  m  is  the  mass  per  unit  area  of  plate. 
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Integrating  Eq.  (4.4)  using  formulas  (4.  1)  and  (4.  3)  yields 


W 


12Mo(X  -  l)/mR2 


(4.5) 


In  Eq.  (4.5),  X  =  p/p  where  p  is  the  pressure  required  to 

8  8 

cause  static  collapse  and  is  [4.2] 


p  =  6M  /R2a2(3-2n) 

8  O 


(4.6) 


a  being  the  ratio  a/R  . 

Now  Eq.  (4.  5)  is  analogous  to  Eq.  (3.  3)  for  the  clamped  beam 
subjected  to  low  pressures.  Thus,  Section  111-2(1)  can  be  taken  over 
to  complete  the  discussion  on  the  simply  supported  plate. 


By  successive  integrations  of  (4.  5),  the  velocity  and  displacement 
of  the  plate  center  are 


W  =  12M  (I  -  I  )/p  mR2 

o  s  r s 

(4.7) 

W  =  12M  (A  -  A  )/p  mR2 

O  B  8 

(4.8) 

where 


I(t)  =  J  p(T)dT  .  Is(t)  =  P|(t 


= j  I(t 


)dT  ,  As  =  pBt  /2 


Motion  ceases  when  W(t)  =  0  ,  at  time  t  =  t£  ,  tay.  Then  from 
(4.  7)  t£  is  determined  by  I(t£)  =  pgt£  .  Consequently  the  maximum 
deflection  is  W(t£)  ,  found  by  substituting  t  =  t£  in  (4.8). 
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*  A,  -  _  ► 


if  1 1 


For  each  a  .  the  peak  load  (xepressnted  nondimcnsionally  as 

\  =  P  / p  )  has  an  upper  bound  X,  such  that  when  it  is  exceeded, 

m  m  *s  1 

the  mechanism  (4.  1)  does  not  apply  because  the  yield  condition  is 
violated.  Determining  X  j  so  that  (4.  3)  is  satisfied  leads  to  the 
following  range  for  Xm  corresponding  to  mechanism  (4.  1) 


lax  *X,  =  2a2(3-2a)/[2oZ(3-2a)-l]  1/2  a  a  a  1 

m  i 

1  a  X  a\,  =  (3  -  2a)/(3  -  4a)  0  s  a  a  1/2 

m  i 


(4.9) 


The  curves  (4.  9)  are  plotted  in  Fig.  4.  4. 


FIG.  4.4  DEFORMATION  MECHANISM  DIAGRAM 


For  peak  pressures  of  blast  pulses  such  that  the  value  of  X  lies 
in  the  range  !  <  X  a  Xj  the  permanent  central  deflection  is  given  by 

W,  =  12M  (A. -A  )/p  mR2  (4.10) 

X  O  X  8  8 

where  W{  =  W(t{)  ,  A{  =  A(tf),  and  t{  =  I(t{)/P8  =  I£/pg  •  This  result 
is  exactly  analogous  to  entire  deformation  by  mechanism  1  for  clamped 
beams  (Section  III). 
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NOMENCLATURE 


radius  of  loading 

r  I(T)dT 

o 

A(t  ) 

Pst2/2 

impulse 

in,) 

p.' 

mass  per  unit  area 

components  of  moment 

fully  plastic  moment 

pressure 

peak  pressure 

static  collapse  pressure 

radial  component 

radius  of  plate 

time 

time  when  motion  ceases 
deflection 
central  deflection 
final  central  deflection 
a/R 

circumferential  coordinate 
components  of  curvature 

p/p, 

an  upper  bound  of  \ 
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SECTION  V 


CLAMPED  RIGID-PLASTIC  PLATES  UNDER  BLAST  LOADING 

A  theoretical  study  is  made  of  clamped  circular  plates  of  rigid- 
plastic  material  subjected  to  blast  loading  uniformly  distributed  over 
the  surface.  The  dependence  of  the  permanent  central  deflection  on 
pressure  and  impulse  is  obtained  when  the  blast  pulse  is  taken  as  a 
rectangular  pulse.  Experiments  are  described  and  the  permanent 
central  deflections  obtained  are  correlated  with  the  theoretical  pre¬ 
dictions. 

1.  Introduction 


The  problem  treated  is  the  response  of  a  clamped  circular  plate 
subjected  to  a  suddenly  applied  pressure  uniformly  distributed  over 
the  whole  of  one  side  of  the  plate.  The  pressure  is  assumed  high 
enough  or  held  on  the  plate  long  enough  to  produce  moderately  large 
plastic  deformations.  Although  the  governing  equations  are  derived 
for  a  general  pressure-time  relationship  they  are  solved  for  the 
simplest  case;  that  of  a  pressure  which  is  held  constant  for  a  time  and 
then  suddenly  released  (a  rectangular  pulse).  It  is  the  variation  of  the 
permanent  central  deflection  with  pressure  and  impulse  (area  under 
pressure-time  curve)  that  constitutes  the  principal  result  of  interest. 

To  simplify  the  analysis  the  plate  material  is  assumed  to  be  rigid- 
perfectly  plas'.ic  in  behavior  obeying  the  Tresca  yield  condition  and  the 
associated  flow  rule.  Only  the  bending  action  of  the  plate  is  taken  into 
account. 

The  response  of  a  clamped  circular  rigid-plastic  plate  to  a 
uniformly  distributed  ideal  impulse  was  found  by  Wang  and  Hopkins  [  5.  l] 
Their  method  of  solution,  established  continuity  and  jump  conditions, 
and  nomenclature  are  employed  here.  Wang  [5.2]  obtained  the  response 
to  an  ideal  impulse  for  simply  supported  plates.  Hopkins  and  Prager. 

[5.  3]  solved  the  problem  similar  to  one  presented  here  but  for  simply 
supported  plates  and  their  results  a/e  used  later  to  assess  the  effect 
of  the  boundary  conditions. 
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A  blast  pulse  may  be  idealized  to  that  shown  in  Fig.  5.  1  with  an 
instantaneous  rise  to  the  peak  pressure  p^  followed  by  a  continuous 
monotonic  decay.  This  property  will  be  used  later.  If  the  peak  pressure 
is  only  slightly  greater  than  the  static  collapse  pressure  p  it  is  reason- 
able  to  expect  the  dynamic  mode  of  collapse  to  be  similar  because  the 
inertia  forces  are  still  small.  In  [5.4]  it  is  shown  that  the  static  collapse 
pressure  is 

2 


p  =  6M  /r 

rs  os 


(5.  1) 


where  Mq  is  the  fully  plastic  moment  per  unit  length  and  rg  is  a 
certain  radius  determined  by  the  equation 


5  +  tn(R/rJ2  *  3(R/rg)2 


(5.2) 

where  R  is  the  radius  of  the  plate.  In  fact,  the  solution  of  (5.2)  is 
rg/R  =  0.73. 

Figure  5.2  shows  the  mechanism  of  deformation  called  mechanism  1, 
corresponding  to  peak  pressures  in  excess  of  p  but  below  a  pressure  p. 
to  be  determined.  The  radius  r^t)  replaces  rg  .  At  the  plate  center 
the  plastic  regime  is  A  in  Fig.  5.3  where  M  =  N  =  Mq  ,  M  and  N 
being  the  radial  and  circumferential  components  of  bending  moment 
(positive  moments  causing  tension  on  the  underside  of  the  plate). 


k  AB  B  BC  C 


«*-«*«•  .it 


FIG.  5.1  IDEALIZED  BLAST  PULSE 


FIG.  5.2  MECHANISM  1 
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At  r  =  r  j (t)  and  r  =  R  the  regimes  are  B  and  C  while  in  the  regions 
0  *  r  *  rj(t)  and  r^t)  *  r  £  R  they  are  AB  and  BC.  A  velocity 
field  which  satisfies  the  flow  rule,  boundary  conditions  and  the  appro¬ 
priate  continuity  and  discontinuity  conditions  [  5.  1  ]  is 

V(l-ar/r  )  0  £  r  £  r  (t) 

wt  =  (5.3) 

Votn  R/r  r  (t)  £  r  £  R 

with 

l/o  =  InR/r  j  +  1  (5t  4) 

In  (5.3),  w  is  the  plate  deflection,  V  is  the  velocity  of  the  plate  center 
and  the  subscript  t  denotes  partial  differentiation.  The  initial  condition 
of  the  plate  at  rest  may  be  expressed  by  V(0)  =  0  . 

The  upper  bound  Pj  of  the  peak  pressure  p^  is  that  pressure 
which  causes  at  the  plate  center  an  inflection  point  in  the  bending 
moment  diagram  for  the  radial  component,  that  is  d2M/dr2  =  0  at 
r  =  0  .  The  inequality  Pm  <  Pj  prevents  violation  of  the  yield  condition 
at  the  plate  center.  This  suggests  that  whenever  the  peak  pressure 
exceeds  p{  a  central  region  0  £  r  £  rQ(t)  <  r^t)  of  the  plate  acquires 
a  uniform  velocity.  This  mechanism,  called  mechanism  Z,  is  shown  in 
Fig.  5.  4  with  the  plastic  regimes  indicated  and  has  the  following  velocity 
field  [5.1] 

V  0  £  r  £  ro(t) 

wt  =  V(l-0(r-ro)/rj)  rQ(t)  £  r  <  r^t)  (5.5) 

VOl n(R/r)  r j(t)  s  r  £  R 


FIG.  5.3  TRESCA  YIELD  HEXAGON  FIG.  5.4  MECHANISM  2 
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’  ^  k*  .* 


where 

l/o  =  tnR/r  j  +  (1  -  rQ/r1)  (5.6) 

By  taking  the  blast  pulse  in  the  form  of  Fig.  5  1  the  value  of 

r  (t)  is  a  maximum  at  time  t  =  0  and  thereafter  decreases  to  zero, 
o'  ' 

During  the  remaining  time  the  plate  deforms  in  accordance  with 
mechanism  1. 

3.  Governing  Equations 

Independent  of  mechanisms  the  equation  of  motion  is 


r 

N  -  (rM)r  s  J  (p  -  mw^Jrdr 


(5.7) 


A  subscript  r  indicates  partial  differentiation  with  respect  to  r  . 

(1)  Mechanism  2 

When  the  peak  pressure  is  large  enough  to  cause  deforma¬ 
tion  by  mechanism  2  (pm  >  Pj)  the  acceleration  to  be  substituted  in 
(5.  7)  is  obtained  by  differentiating  (5.  5)  with  respect  to  time.  The 
circumferential  component  of  the  bending  moment  N  is  eliminated 
by  using  the  yield  condition  (Fig.  5.  3).  Due  to  the  three  properties 
M  =  N  =  0  in  0  <  r  *  rQ  ,  M(rltt)  =  0  ,  and  M(R,  t)  =  -  Mq  , 
carrying  out  the  integration  in  (5.  7)  provides  the  following  three 
equations. 

V'  =  (5.8) 

ry 

V#(?+T})n[2?(3-3T]  +  ii2)  +  ti(6^8ti  +3ti2)]  -  V?  rft  ?(6-8ti+3ti2)+  ti(1-ti)(4-3t|)] 

-VVn2[25(3-2n)  +  n(4-3n)]  =  (3-3  n  +  ri2)  -  1]  e2^(?  +  Ti)2 

(5.9) 

V,(§  +  r1)[3e2?-3  -2?(3-3ri+3v|2-r?)]-V§/[3e2?-3  -2?|3-r2(l-n)(3-2n)} 

-  2»2(3-6n+6n2-2r3)l-Vn#[3e2?-3  -2§  (3  -  3t}2+ 2ti3)  -  6§2(1-t/j] 

=  [3\e2(?i  ?)(e25-l)/2-(l+?)]e2^  +  n)2 

(5.10) 


* 


r  jflGr 
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The  new  variables  that  have  been  introduced  in  the  derivation  of  (5.8), 
(5.  9)  and  (5.  10)  are  defined  by  [  5.  l] 

?  =  {-nfl/pj)  T|  —  1  —  Pq/Pj  Pj  =  fj/R 

Pc =  r0/R  X  =  p/pg  ?g=tn(l/ps)  pg  =  rg/R 

Also,  the  primes  indicate  differentiation  with  respect  to  the  variable 
T  '  where 

T  =  12M  t/mR2 
o 


The  value  of  r  ,  and  hence  p  and  §  ,  is  the  solution  of  Eq.  (5.2) 

8  8  8 

while  the  static  collapse  pressure  pg  is  given  by  (5.  1). 

(2)  Mechanism  1 


When  the  value  of  X  =  p/pg  falls  in  the  range  1  £  X  s  Xj 
where  Xj  =  Pj/pg  deformation  occurs  by  means  of  mechanism  1.  The 
governing  equations  can  be  found  by  substituting  in  (5.  7)  the  accelera¬ 
tions  found  by  differentiating  (5.  3)  with  respect  to  time  and  carrying 
out  the  integration  noting  that  M(r^,t)  =  0  and  M(R,  t)  *  -  . 

Alternatively,  they  can  be  found  by  setting  r)  =  l(pQ  =  0)  and  p  '  =  0 
in  Eqs.  (5.  9)  and  (5.  10).  Performing  either  of  these  operations  yields 

V'(5+  1)(2?+  1)  -  VS'5  =  [xe  8  -  l]e2?(?+l)2  (5.11) 


V'(?+l)(3e2?-3-45)-V§#(3e2?-3-6?-2^)  =  (3Xe  8  (e2?-l)/2 

-  (l+$)]e2?(5  +  l)2 


2Cg 

Whenever  V  '  =  Xe  / 2  an  inflection  point  in  the  bending  moment 

diagram  for  M  occurs  at  the  plate  center  which,  for  pulses  of  the 

type  shown  in  Fig.  5.  1,  occurs  immediately.  This  condition  gives 

X,  where  X .  =  P  /p  . 

1  1  m  *s 


(5.12) 
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(1)  Mechanism  2  Xj  <  \ 

Specialising  to  a  rectangular  pulse  a  solution  of  Eqs. 

(5.8),  (5.9)  and  (5.  10)  is  obtained  if  it  is  aiisumed  that  §  and  r]  are 
constants  while  the  load  is  acting.  This  means  that  while  p  is  constant, 
rQ  and  r^  are  constants.  This  phase  of  motion  will  be  called  phase  la. 
Setting  5  '  -  V  =  0  in  (5.  9)  and  (5.  10),  and  substituting  V  '  from  (5.  0) 
in  (5.  9)  and  (5.  10)  gives 

2(5.-?)  3 

2(5+0  =  Xe  riP-n)  (5.  13) 

2(5  *5)  Zl  Z  3 

2(?+r1)(l  +  S)  =  Xe  1  [3e^(5-l  +  T!)  +  §(3-6^+611-2x0  +  3(1-^)] 

(5.  14) 

The  lower  bound  Xj  of  X  may  be  found  by  substituting 
n  =  1  ( p  *  0)  in  (5.  13)  and  (5.  14).  Doing  this  gives 

2?.  2? 

Xje  s=2(?+l)e^  (5.15) 

where  5  determined  by  the  equation 

3?e25  =  1  (5.  J 6) 

From  (5.  15)  and  (5.16),  X  j  ~  1.  998  and  5  =  0.2163  (p^  =  0.805). 

For  a  given  value  of  X  such  that  X  >  Xj  Eqs.  (5.  13)  and 
(5.  14)  give  the  initial  values  of  5  and  *1  •  Some  numerical  valuer 
are  shown  in  Table  5.  1. 

Let  the  pulse  end  at  time  t  =  t  or  when  t  '  -  t  ' .  If 

o  o 

the  velocity  of  the  plate  center  at  this  time  be  denoted  by  Vy  ,  integra¬ 
tion  of  Eq.  (5.  8)  gives 

2? 

V  =  1/2  Xe  *t  '  =  I/m 
o  o 

where  I  =  pt  is  the  impulse, 
o 
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Table  5.  1 

INITIAL  VALUES  OF  %  AND 


Mechanism 


1.  1 
1.2 
1.  3 
1.4 
1.  6 
1.8 
1-9  „ 
1.  998 

2.0 
2.24 
3.80 
6.  51 
9.48 
15.71 
32.  53 
56.  13 


.  302 
.290 
.279 
.268 
.248 
.232 
.224 
.2163 

.216 
.2 
.  14 
.  10 
.08 
.  06 
.04 
.03 


(0) 

V 

p/°> 

•o'0’ 

1.0 

.  739 

0 

1.  0 

.  748 

0 

1.0 

.  756 

0 

1.  0 

.  764 

0 

1.  0 

.  780 

0 

1.  0 

.793 

0 

1.0 

.  800 

0 

1.0 

.805 

0 

.998 

.806 

.  001 

.883 

.819 

.  096 

.569 

.869 

.375 

.395 

.905 

.547 

.313 

.923 

.  634 

.233 

.  942 

.722 

.  156 

.  961 

.811 

.  116 

.970 

.858 

Value  of  \ 


Now  V  =  pt/m  ,  so  by  integration,  the  central  deflection 


o{  the  plate  at  time  t  =  t  is 

o 


2  7  7  2§ 

6o  =  1  ^2mP  =  (lR/  12mMQ)/Xe  8 


(5.17) 


When  t  >  tQ  ,  no  pressure  acts  on  the  plate  so  that  X  =  0 
and  hence  from  Eq.  (5.8)  V'  =  0  .  This  means  that  the  central  region 
of  the  plate,  0  s  r  s  rQ(t)  moves  at  a  constant  velocity.  It  is  evident 
from  (5.9)  and  (5.  10)  that  %  and  n  can  no  longer  be  treated  as  constants. 
Introducing  a  dimensionless  time  t  =  12M  (t-t  )/mR2V  =  12M  (t-t  )/lR2 

O  O  o  o  *  o* 


Eqs.  (5.  9)  and  (5.  10)  become 

?  '  t  ?(6-8r|  +  3rj^)  +  ri(l-ri)(4-3ri )]  +  n '  [  2?(3 -2r| )  +  ri  (4-3n )]  =  e2  1?+r)2/r  2 

(5.  18) 

?/[3e2?-3-25{3-2T12(l-T1)(3-2T1)}  -2?2(3-6ti+  6«n2-2n3)] 

1  J  (5.  19) 

+  T1'l3e25-3-2$(3-3Ti2+2T13)  -  6?2(l-n)2]  =  e2^+ n )2  (§+ 1) 

where  now  the  primes  denote  differentiation  with  respect  to  i.  In 
(5.  8),  (5.  9)  and  (5.  10)  derivatives  are  with  respect  to  t  '  (dT  '  =  V  dT  ). 

This  phase  of  the  motion  described  by  Eqs.  (5.  18)  and  (5.  19),  called 
here  phase  lb,  is  analogous  to  the  first  phase  of  motion  in  the  impulse 
problem  [  5.  l],  except  that  the  initial  values  of  ?  and  r\  are  ?^  and 
rf°^  from  Eqs.  (5.  13)  and  (5.  14)  instead  of  ?  -  1  and  r\  =  0  .  Phase  lb 
ends  when  r\  *  1(po  *  0). 

An  outline  of  the  method  used  to  solve  Eqs.  (5.  18)  and  (5.  19), 
which  is  essentially  the  same  as  that  used  in  [  5.  l],  follows  immediately. 

F  rom  (5.  18)  and  (5.  19) 

d?/dn  =  -P(?»  ri  )/Q(5  *  +i  ) 

where 

P(5,ti)  =  3(1+2? +2?2-e2?)  +  T1U2(4-3T1)(?+l)  -  4?2(ti2-3*  +  3)] 

Q(?,-n)  =  3(1+2?+  2?2-e2')  +  r\[  t12(4-3t1  )(l-n)  +  ?  t12(6-8t1+3t12)-3?2(4-6-  +4^2-r3)] 

Let  d?/dri  =  -P(§^°\  T|^)  /Q(?^°\  =  mQ  ,  where  ?^°*  and  are 

determined  from  (5.  13)  and  (5.  14)  and  depend  only  on  \  .  Choose  a  new 
value  of  ?,  say  ?^  slightly  greater  than  ?^  ;  this  corresponds  to  the 
assumption  that  decreases  in  phase  lb.  Then  the  new  value  of  r,  , 
say  •  *8  determined  by  v ^  +  (?^  -  ?^)/mQ  •  Now  find  a 

new  value  for  the  ratio  d?/dn  from  -P(?^,  n^)/Q(?^i  t/1^)  *  • 

Continue  this  process  until  t]  =  1  ,  when  the  plastic  hinge  circle  reaches 
the  center  and  phase  lb  terminates.  The  remaining  deformation  occurs 
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in  mechanism  1.  Figure  5.  5  shows  the  ?  -  i|  trajectories  obtained  by 
this  procedure  for  eight  values  of  \  .  To  find  the  time  at  the  end  of 
phase  lb  first  find  the  function  §'(§  ,  r, )  using  (5.  18)  and  (5.  19).  Then 
the  numerical  value  of  %  '  ^  is  ?  '  ,  r>^).  Now  find  the  time 

increment  from  t  ^  =  (5^  -  ?^)/?/^.  Continue  this  process 

along  a  §  -  r)  trajectory  until  r|  =  1  .  The  sum  of  the  time  increments 
then  gives  "  ^  corresponding  to  the  duration  of  phase  lb. 


FIG.  5.5  VARIATION  OF  £  WITH  rj  DURING  PHASE  1b 


If  6  j  is  the  deflection  of  the  plate  center  at  the  end  of 
phase  1  it  is  obtained  from  the  equation 

\  -  &0  =  Vo(tl  '  V  =  ^V1 2mMo  <5-20> 

in  which  6Q  is  given  by  (5.  17). 

The  final  phase  of  the  motion,  called  phase  2,  is  governed 
by  Eqs.  (5.9)  and  (5.  10)  with  tj  =  1  .  tj  '  =  0  and  X  =  0  .  With  these 
special  values  and  retention  of  the  dimensionless  time  variable  t  used 
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in  phase  lb  Eqs.  (5.  9)  and  (5.  10)  become 

C'(5+  U(2?+  1)  "  CS'S  =  -(?+l)2«2?  (5-21> 

£'(?  +  l)(3eZ?  -3-4?)  -  £§  '(3e2?  -3-6?-2?2)  =  -('  +  l)’e2?  (5.  22) 

where  f=V  V  .  Equations  (5.  2  1)  and  (5.  22)  can  be  derived  directly 
3  o 

by  substituting  the  velocity  field  defined  by  (5.  3)  and  (5.  4)  in  the  equation 
of  motion  (5.  7). 


From  (5.21)  and  (5.22) 


c  =<^,exp[-/ 77^^17?]  (5‘ 


23) 


in  which  is  the  value  of  ?  at  the  end  of  phase  lb. 

Motion  ceases  when  £  =  0  and  this  occurs  when  §  =  ~  478 

which  is  the  solution  of  the  equation 

4  +  7?  +  2?2  -  3eZ?  =  0 

Figure  5.  6  shows  the  variation  of  £  with  ?  . 

Let  t2  be  the  value  of  t  when  motion  ceases.  Then 


?- 


t2 


?2  13 2 

-  t,  =  f  4  =  f  (3e^-4^  -6|  -3)£d| 

1  J  5  J  e2?(?+  l)(4  +  7£+2?  -3e25) 


(5.24) 


Then 


2?  _ 4?2  - 

«i  ?i 

Finally,  let  the  central  deflection  when  t=  T2(t=t2)  be  62  . 

T  •, 


A  A  -  IR 

62  “  61  "  TZm, 


/vd,  =  f 

J  12mMQ  y  eZ^(?+l)(4+7?  +  2?Z-3e 

(5.25) 


“ 
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/) 


The  value  of  is  almost  constant  during  phase  2.  As  \  increases 
the  diagrams  become  progressively  stronger  in  resemblance  to  that 
for  an  ideal  impulse  (see  Fig.  3  of  [  5.  l] ). 


FIG.  5.9  VARIATION  OF  p0  AND  p]  WITH  r.  A  =  56.13 


i 
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(2)  Mechanism  1  1  <  X  £  Xj 

When  the  constant  pressure  lies  between  pg  and  p^  the 
equations  governing  the  motion  during  phase  la  can  be  obtained  by 
substituting  in  (5.9)  and  (5.  10)  the  special  values  r|  =  1  ,  r  '  -  0  and 
§/  =  0  .  These  substitutions  give 


,  r  2(5  -§)  ,  2? 

V'  (2§+  1)  =  [xe  s  -  l]e^(?+  1] 


(5.26] 


2(S  -5) 


(3e2?-3 -4§)  =  (3Xe  8  "  (e2?- l)/2  -  (•  +  l)]e2?(?  +  1)  (5.27) 


where  the  prime  denotes  differentiation  with  respect  to  t  '  .  Equations 
(5.  26)  and  (5.  27)  also  result  when  the  velocity  field  defined  by  (5.  3) 
and  (5.4)  is  substituted  in  the  equation  of  motion  (5.  7)  and  use  is  made 
of  the  properties  M(rj,t)  =  0  and  M(R,  t)  =  -  Mq  . 

Combining  Eqs.  (5.26)  and  (5.27)  to  eliminate  V;  provides 
the  following  equation  which  determines  ?  once  X  is  given  (by  choosing 
5's  and  finding  X's  a  X  -  ?  curve  can  easily  be  obtained): 

r  2(?s_?)  25  ,  r  ,  25 

[3Xe  8  (e^  -  l)/2  -  (§  +  1)](2?+1)  =  [  Xe  8  l](3e^-3-4§) 


(5.28) 


Knowing  §  from  (5.28)  the  velocity  of  the  center  of  the 
plate  by  integration  of  (5.26)  is 


r  2(f  -?) 

V  -  l  Xe  -  l]e^  (?+  1)t7(25+  1] 


(5.29) 


The  central  deflection  of  the  plate  at  time  t  =  t  is 


(2?+  l)Xe 


Phase  2  is  essentially  the  same  as  the  phase  2  described 
earlier  when  it  followed  a  phase  1  deformation  according  to  mecha¬ 
nism  2. 


From  (5.  29)  , 

V. 


-  21 


£1  [1  - 
m 


"ZTT 


Xe 


T2  ”  T1  ^w^ere  now 


Now  set  C=  V/V  as  was  done  earlier.  Then  C  and 

°  2? 

=  2/Xe  *  because  there  is  no  phase  lb) 


T1  =  To 


are  given  by  (5.23)  and  (5.24),  and  in  place  of  (5.25),  the  central 
deflection 

r.  c 

>2 


6-  when  motion  ceases  is  given  by 

2  §5 


2  2 

ft  -  IR 
°2  1  ‘  imM 

o 


[  i-  -nfrrrl2  &  f 

XeZ"s  Z?  +  1  *  e  ?(?+l)(4+7?+2f  -3e  5 


) 


(5.31) 


In  (5.31),  6j  =  6q  (no  phase  lb)  which  is  given  by  formula  (5.30).  Also, 
is  the  solution  of  Eq.  (5.28). 


5.  Conclusions 

Figure  5.  10  shows  the  relationship  among  the  final  central 
deflection  6,  the  impulse  per  unit  area  I  and  the  pressure  pm  in 
nondimens ional  form  for  both  the  clamped  and  simply  supported 
plates.  The  curve  for  the  clamped  plates  is  obtainable  from  formulas 
(5.17),  (5.20)  and  (5.25)  fo^  X  >  Xj  and  from  formulas  (5.  30)  and 

(5.  31 »  for  1  <  X  s  Xj  .  That  for  the  simply  supported  plates  is  obtain¬ 
able  from  the  results  in  [5.3].  In  using  Fig.  5.  10  it  should  be  noted 
that  p  =  6M  It  for  the  clamped  plates  and  p  =  6M  /R  for  the 

‘8  OS  S 

simply  supported  plates.  The  static  collapse  pressure  for  the  former 
is  1.875  times  tht.t  for  the  latter. 
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FIG.  5.10  RELATIONSHIP  BETWEEN  8,  A  AND  I 

Figure  5.  11  shows  the  relationship  that  must  exist  between  the 
pressure  and  impulse  to  provide  a  given  central  deflection  of  a  clamped 
plate.  (The  corresponding  curve  for  a  simply  supported  plate  lies 
almost  on  top  of  the  shown  curve. )  In  other  words,  points  on  the  curve 
define  a  family  of  rectangular  pulses  each  member  of  which  produces 
the  same  central  deflection  of  a  clamped  plate.  The  coordinates  of 
Fig.  5.  11  have  been  rendered  nondimensional  by  using  \  and  Ij  , 
the  latter  being  the  ideal  impulse  which  causes  the  same  central  de¬ 
flection  as  the  rectangular  pulses.  In  fact,  6  =  0.  07  I?R2/mMc  [5.  l], 
a^d  6  =  vl  R  /mM  where  v  is  obtainable  from  Fig.  5.  10  (v  depends 


FIG.  5.11  PRESSURE-IMPULSE  DIAGRAM 


From  Figs.  5.  10  and  5.  11  the  following  conclusions  can  be  drawn: 

1.  For  a  given  impulse,  the  central  deflection  6  increases 
with  the  pressure  pm  becoming  a  maximum  when  the 
pressure  is  infinite, that  is,  when  the  given  impvise  is 
applied  as  an  ideal  impulse.  This  can  be  seen  in  Fig.  5.  10 
for  both  clamped  and  simply  supported  plates  by  fixing  I  , 
varying  X  and  observing  v  . 

2.  Again  for  a  given  impulse,  consideration  of  how  6  increases 
with  X  in  Fig.  5.10  reveals  that  for  clamped  and  simply 
supported  plates  respectively  ever  85%  and  90%  of  the  maxi¬ 
mum  central  deflection  (6  for  impulse  Ij)  are  obtained  when¬ 
ever  the  pressure  is  greater  than  six  times  the  static 
collapse  pressure,  that  is,  whenever  X  >  6  .  Below  this 
value  of  X  the  decrease  of  6  with  X  is  quite  pronounced. 

3.  For  a  given  central  deflection  6  ,  Fig.  5.  11  shows  that  as 
the  pressure  is  decreased  from  infinity  to  a  value  corres¬ 
ponding  to  about  X  =  6  the  increase  in  impulse,  over  the 
ideal  impulse,  neceasary  to  maintain  the  given  deflection 
6  is  less  than  7%.  (1  <  I/Ij  <  1.  07  for  6  <  X  <  CO  )  A 
large  increase  in  impulse  is  required  to  maintain  6  as 

X  decreases  further,  especially  in  the  range  1  <  X  <  2  . 
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4.  In  comparing  the  effect  of  the  boundary  conditions  the 

static  collapse  pressure  for  clamped  and  simply  supported 
edges  are  in  the  ratio  1.875  :  1  .  From  the  results  of 
[5.  1.  5.2]  the  ratio  of  ideal  impulses  to  achieve  the  same 
central  deflection  is  1.34:  1  .  This  latter  ratio  is  the 
minimum  value  for  rectangular  pulses  for  all  values  of 
\>1  which  implies  considerable  strengthening  due  to 
clamping  the  edges. 

6.  Description  of  Experiments 


Circular  plates  of  6061-T651  aluminum  and  of  various  thicknesses 
ranging  from  3/16  inch  to  3/8  inch,  clamped  at  either  8  inches  or  12 
inches  diameter,  were  subjected  to  the  same  type  of  pulse  generated 
in  an  oxyacetylene  shock  tube  32  inches  long  and  having  the  same 
internal  diameter  as  the  supports.  The  experimental  arrangement  is 
shown  assembled  in  Fig.  5.  12  and  dismantled  in  Fig.  5.  13.  Thick 
steel  annular  supports  provided  clamping  against  rotation  only,  a 
spacer  keeping  their  distance  apart  slightly  greater  than  the  plate 
thickness.  Apart  from  circumferential  membrane  forces  and  unavoid¬ 
able  frictional  effects  the  plate  material  was  allowed  to  displace  freely 
inwards  during  deformation.  To  ensure  good  clamping  the  overall 
diameters  of  the  plates  tested  were  9  1/2  and  15  inches  but  to  minimize 
the  build-up  of  membrane  forces  they  were  slotted  radilly  at  1/2-inch 
intervals  around  the  rim.  The  length  of  the  i  lots  were  such  that  their 
ends  did  not  pass  over  the  support  circle  during  deformation.  Gas 
seals  of  heavy  silicon  grease  were  provided  at  the  junctions  of  the  upper 
steel  support  with  both  the  plate  and  the  shock  tube  to  reduce  to  rare¬ 
faction  of  the  applied  shock  wave  due  to  relief  of  pressure  at  the  junctions. 

The  pulse,  shown  in  Fig.  2.4,  was  taken  to  be  that  recorded  under 
a  similar  configuration  by  a  rod  gage  and  is  therefore  the  pressure 
acting  on  a  fixed  target.  This  implies  the  assumption  that  the  modifi¬ 
cation  of  the  pulse  due  to  plate  motion  is  small.  The  peak  pressure 
Pm  is  obtained  from  the  pressure  record  of  Fig.  2.  4  by  extrapolating 
smoothly  the  mean  pressure  curve  back  to  the  starting  time.  This 
results  in  a  peak  pressure  of  pm  =  1440  lb/in  .  Now  the  mean  pressure 
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.  v.c, 


FIG.  5.12  EXPERIMENTAL  ARRANGEMENT 


curve  is  fitted  by  the  equation  p  =  pme”8t  to  find  s  .  By  integration, 

I  =  p  /•  and  hence  the  impulse  (area  under  pressure -time  curve)  is 

m  2 
known.  The  value  of  the  impulse  is  I  =  0.335  lb-sec/in  . 


■-i. 


FIG.  5.13  EXPERIMENTAL  ARRANGEMENT  (DISMANTLED) 


Tensile  tests  performed  with  the  Instron  machine  on  standard 
ASTM  specimens  taken  with  and  across  the  direction  of  rolling  of  the 
plate  material  provided  stress-strain  diagrams  from  which  the  yield 
stress  was  determined.  The  yield  stress  is  here  defined  as  the  stress 
at  the  point  of  intersection  of  the  straight  line  approximations  to  the 
elastic  and  strain-hardening  portions  of  the  stress -strain  diagram.  An 
average  of  the  yield  stresses  obtained  from  specimens  with  and  across 
the  rolling  direction  was  taken.  The  plate  material,  6061-T651 
aluminum,  was  chosen  because  it  exhibits  only  slight  strain-hardening 
and  is  insensitive  to  strain  rate. 

After  the  blast  loading  tests  the  permanent  central  deflections 

6  were  measured, 
ex 
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7.  Experiment*!  Results  and  Observations 

The  results  of  the  experiments  described  above  are  presented 
in  Table  5.2  and  Fig.  5.  14.  Correlation  with  the  predictions  of  the 
rigid-plastic  theory  takes  place  through  the  ratio  &ex/6th  of  the  central 
deflections  in  Table  5.2  and  by  means  of  the  X  versus  IR/(12mMQ6) 
plot  of  Fig.  5.  14  which  is  a  form  of  p-I  diagram. 

The  main  observation  is  that  the  central  deflection  ratios  6ex/fith  • 
by  ranging  from  0. 14  to  0,  38,  show  that  the  theory  overestimates  the 
central  deflections  from  the  present  experiments  by  factors  ranging 
from  about  3  to  7.  That  the  plates  are  much  stronger  than  predicted 
is  due  to  the  following  reasons: 

(1)  In  the  rigid-plastic  theory  the  effects  of  elastic  strain 
energy  and  vibrations,  strain -hardening,  and  strain  rate 
are  neglected  all  of  which  add  strength  to  a  structure.  How¬ 
ever,  the  material  used  has  little  strain -hardening  and  is 
fairly  insensitive  to  strain  rate. 

(2)  If  the  deflections  become  large  enough  membrane  forces 
become  significant  and  the  theory  considers  a  bending  action 
only.  In  general,  it  can  be  seen  in  Table  5.2  that  as  the 
ratio  6ex/R  increases  the  ratio  6ex/6tjj  decreases. 

(3)  Frictional  and  circumferential  membrane  forces  near  the 
plate  rim,  or  rather  near  the  circle  passing  through  the 
ends  of  the  radial  slots,  add  strength  to  the  plate. 

(4)  The  experimental  pulse  was  exponential  whereas  the 
theoretical  pulso  was  rectangular  and  for  a  given  peak 
pressure  and  impulse  the  latter  predicts  larger  deformations. 

(5)  The  applied  impulse  is  probably  less  than  that  recorded 
against  a  rigid  target.  In  Fig.  5.  14  all  the  experimental 
points  would  be  moved  to  the  left  (in  the  direction  of 
decreasing  impulse)  and  hence  into  closer  agreement  with 
the  theoretical  curve. 


Table  5.2 

EXPERIMENTAL  RESULTS 


Thickness 

d 

(inch) 

Yield 

Stress 

0 

(lb /in2) 

Static 

Collapse 

Pressure 

P.2 
(lb /in'  ) 

k 

P. 

3/8 

43.  350 

254 

5.67 

5/16 

42,700 

174 

8.  31 

1/4 

41.  800 

109 

13.22 

3/16 

39.  300 

58 

25.05 

3/8 

44.  550 

688 

2.  10 

5/16 

42,700 

391 

3.68 

1/4 

42,200 

248 

5.81 

3/16 

39,  300 

129 

11.10 

Peak  preaaure 
Impulse 
Mata  density 
Matt 

Plastic  moment 


1440  lb /In 

0.  335  lb-aec/ln2 

0.  OOOZ53  lb-tec2/in4 

pd  lb-sec2/in* 

o  d2/4  lb-ln/in 
o 


Central 


ex 

(inches) 


1R 

‘ex 

nr 

J  lsirr  •  fc 

y  o  ex 

2.72 

052 

2.61 

.  057 

2.  61 

.  057 

2.29 

.  13 

2.44 

.  12 

2.59 

.20 

2.56 

.21 

2.66 

.  19 

2.  56 

.31 

2.56 

.31 

2.32 

.047 

2.  i 5 

.055 

2.  I  7 

.094 

2.  17 

.094 

2.30 

.17 

2.  30 

.  17 

2.  90 

.26 

Material:  6161-T651  aluminum 
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FIG.  5.14  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


The  duration  of  the  fixed  target  pulse  is  about  800  psec  (Fig.  2.  3) 
and  much  longer  than  the  fundamental  response  times  of  quarter 
periods  T/4  of  the  fundamental  modes  of  vibration  of  the  plates  (see 
Table  5.  3  below).  It  is  known  from  the  photodiode  experiments 
(Section  II)  that  both  the  shape  and  impulse  of  the  rigid-target  pulse 
is  affected  by  the  target  mobility.  An  additional  complication  is  that 
the  target  mobility  varies  over  the  plate.  This  interaction  is  important 
and  requires  further  study. 
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Table  5.  3 


RESPONSE  TIMES  OF  FUNDAMENTAL  MODES 


Radius 

Thickness 

Response  Time 

a 

d 

T/4 

(inches) 

(inch) 

(psec) 

NOMENCLATURE 


plate  thickness 

impulse 

ideal  impulse 

mass  per  unit  area 

components  of  bending  moment 

fully  plastic  moment 

pressure 

mechanism  pressure  bound 
peak  pressure 
static  collapse  pressure 
radial  coordinate 

plastic  radii  of  regime  boundaries 

plate  radius 

time 

pulse  duration 

time  at  end  of  phase  lb 

time  at  end  of  phase  2 

fundamental  period  of  plate  vibration 

velocity  of  plate  center 

v(to) 

plate  deflection 

plate  central  deflection  at  time  tQ 
plate  central  deflection  at  time  t^ 
plate  central  deflection  at  time  t ^ 
experimental  final  central  deflection 
theoretical  final  central  deflection 
V/V 

i  -  p„/p, 

p/p, 

Pl'P.  ‘  Pm/p. 
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NOMENCLATURE  (Continued) 


? 

5. 

Si 

^2 


T1 

T2 


/ 

T 


/ 


0 


6Mm/lV 

o 

tntl/p^ 


^n(l/pg) 

§  at  time  tj 
5  at  time  t^ 


VR 

r,/R 

ra/R 

12M0(t-t0)/IR2^ 

12M  (t.-t  )/lRfc 

12M°(t2-t  )/IR2 

12Mot/mR“ 

12M  t  /mR2 
o  o  . 

(inR/rj  +  if1 


► 
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SECTION  VI 


CLAMTJED  CIRCULAR  RIGID-PLASTIC  PLA  f ES 
UNDER  BLAST  LOADING  (Partial  Loading) 

A  theoretical  study  is  made  of  clamped  circular  plates  of  rigid- 
plastic  material  subjected  to  blast  loading  uniformly  distributed  over 
a  central  circular  area.  The  dependence  of  the  permanent  central 
deflection  on  pressure,  impulse  and  loaded  area  is  obtained  when  the 
blast  pulse  is  taken  as  a  rectangular  pulse.  Experiments  are  described 
and  the  permanent  central  deflections  obtained  are  correlated  with  the 
theoretical  predictions. 

1.  Introduction 


The  problem  treated  is  the  response  of  a  clamped  circular  plate 
subjected  to  a  suddenly  applied  pressure  uniformly  distributed  over  a 
central  circular  area.  The  pressure  is  assumed  high  enough  or  held  on 
the  plate  long  enough  to  produce  moderately  large  plastic  deformations. 
Although  the  governing  equations  are  derived  for  a  general  pressure¬ 
time  relationship, they  are  solved  only  for  the  simplest  case;  that  of  a 
pressure  held  constant  for  a  time  and  then  suddenly  released  (a  rec¬ 
tangular  pulse).  It  is  the  variation  of  the  permanent  certral  deflection 
with  pressure,  impulse  (area  under  the  pressure -time  curve)  and  the 
loaded  area  that  constitutes  the  principal  result.  Because  of  interest 
in  moderately  large  plastic  deformations  and  in  order  to  simplify  the 
analysis  the  plate  material  is  assumed  to  be  rigid-perfectly  plastic 
obeying  the  Tresca  yield  condition  and  the  associated  flow  rule.  Only 
the  bending  action  of  the  plate  is  taken  into  account. 

Past  work  on  the  dynamic  response  of  a  rigid-plastic  circular 
pl^te  concerns  blast  or  impulsive  loading  uniformly  distributed  over 
the  entire  area.  In  Section  V  a  clamped  plate  subjected  to  a  rectangular 
pulse  is  treated.  Wang  an'l  Hopkins  [5.  l)  found  the  response  of  a 
clamped  plate  to  an  ideal  impulse.  A  similar  method  of  solution,  the 
established  continuity  and  jump  conditions,  and  similar  nomenclature 
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are  employed  here  as  is  done  in  Section  V.  The  responses  of  a  simply 
supported  plate  to  an  ideal  impulse  and  a  rectangular  pulse  were  found 
by  Wang  [5.2)  and  by  Hopkins  and  Prager  [  5.  3]  respectively. 

The  present  problem  is  quite  similar  to  that  of  Section  V  but, 
due  to  the  loading  covering  only  a  central  part  of  the  plate  (Fig.  6.  1), 
two  additional  mechanisms  of  deformation  can  exist.  An  analogous 
situation  exists  for  beams  and  the  two  additional  mechanisms  are  2b 
and  3  of  Section  ITT. 

2.  Mechanisms  of  Deformation 

A  blast  pulse  may  be  idealized  to  the  form  shown  in  Fig.  6.2 
with  an  instantaneous  rise  to  the  peak  pressure  p^  followed  by  a 
continuous  monotonic  dec?.y.  It  is  assumed  that  this  form  of  pulse 
allows  a  steady  progress  through  the  various  mechanisms  of  deformation 
that  v  4llbe  described. 


p(t) 

i  ■> 

j •  I 

FIG.  6.1  CIRCULAR  PLATE  PROBLEM 


FIG.  6.2  IDEALIZED  BLAST  PLATE 


During  deformation  the  plate  is  divided  into  annular  regions  in 
each  of  which  exists  a  certain  plastic  regime  defined  by  a  vertex  or 
side  of  the  Tresca  yield  hexagon  (Fig.  6.3)  relating  M  and  N  ,  the 
radial  and  circumferential  components  of  bending  moment  (positive 
moments  cause  tension  on  the  underside  of  the  plate),  in  this  problem 
only  the  portion  ABC  is  involved.  Associated  with  these  annular  regions 
are  velocity  fields  which  must  satisfy  the  flow  rule,  boundary  conditions, 
and  the  appropriate  continuity  and  discontinuity  conditions  [5.  1 ) . 
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may  be  expressed  in  the  form  [  5 


From  the  equation  of  equilib¬ 
rium,  the  static  collapse  pressure 
pg  ,  acting  on  a  central  circular 
area  of  radius  r  =  a  of  a  plate  of 
radius  r  =  R  ,  is  found  by  assuming 
that  in  the  circle  0  £  r  *  rg  ,  the 
plastic  regime  is  AB,  A  being  at 
the  plate  center  and  B  at  radius 
r  =  rg  ,  while  in  the  outer  annular 
region  r  £  r  £  R  the  plastic 
regime  is  BC  ,  C  being  at  the 
plate  support  r  =  R.  The  result 
4] 


PjjR2a2(l -2a/3og)/2MQ  =  1  0  <  a  *  pg  (6.1) 


where  pg  is  the  solution  of 


2a(l +6n  l/pg)/3ps  =  1  (6.2) 

and 

PbR2P2/6Mq  =1  Ps  s  a  s  1  (6.3) 

\€ 

where  now  p  is  the  solution  of 

s 

p2(5+24nl/pg)  =  3a2(l  +  2tn  1/a)  (6.4) 

In  Eqs.  (6.1),  (6.2),  (6.  3)  and  (6.  4),  p  =  r  /R  ,  a=a/R  and  M 

ss  o 

is  the  fully  plastic  moment. 

For  values  of  the  peak  pressure  in  excess  of  the  static  collapse 
pressure  the  plate  deforms  initially  in  one  of  four  mechanisms  depending 
on  the  parameters  \  *  p  / p  and  a  =  a/R  . 
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(a)  Mechanism  3 


For  large  value*  of  \  it  is  assumed  that  the  initial  mech¬ 
anism  of  deformation  consists  of  a  central  portion  of  radius  rQ(t) 
moving  at  a  velocity  V(t)  ,  an  outer  annular  portion  r2(t)  «  r  *  R  at 
rest,  and  an  annular  region  between  undergoing  plastic  deformation 
(Fig.  6.4).  At  r  =  rQ(t)  and  r  =  r2(t)  (rQ  <  r2)  two  plastic  hinge 

circles  exist  and,  for  blast  pultes 
of  the  kind  shown  in  Fig.  6.  2,  the 
inner  radius  is  assumed  to  decrease 
until  the  center  is  reached  while  the 
outer  increases  until  the  support  is 


O'  _  _* 

*  r* _ L 


».  ». 


Uf  A  A 


I  I  I 


I  J|il 

I  A  | AB | BC I C I 

!  till 


U-4WI-H 


reached.  At  the  radius  r  =  r^t) 
the  plastic  regime  B  exists  and  the 
rate  of  the  radial  component  of 


FIG.  6.4  MECHANISM  3  curvature  changes  sign  across  this 

circle.  The  plastic  regimes  existing 
elsewhere  are  shown  in  Fig.  6.  4. 

The  velocity  field  meeting  the  above  description,  satisfying  the  flow 
rule,  boundary  conditions  and  the  appropriate  continuity  and  disconti¬ 
nuity  conditions  [5.  l]  is 


w. 


■ 


V  (  1  -0(r  -  rQ)/r1) 


V  otn  r2/r 


0  s  r  s  rQ(t) 
r o<t)  s  r  s  r ^(t) 
i* i (t)  *  r  «  r2(t) 
r2(t)  «  r  s  R 


(6.5) 


where  l/o  =  Inrj/rj  f  1  -  rQ/rl  . 

In  (6.  5),  w  is  the  plate  deflection  and  the  subscript  t 
denotes  partial  differentiation. 

Deformation  proceeds  under  another  mechanism  once  one 
of  the  hinge  circles  reaches  its  terminal  position.  If  the  outer  hinge 
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circle  reaches  the  support  first  the  ensuing  mechanism,  called  mech¬ 
anism  2a,  is  described  by  the  velocity  field  (6.  5)  with  r^  =  R  •  If  the 
inner  hinge  circle  reaches  the  center  first  the  ensuing  mechanism, 
called  mechanism  2b,  is  described  by  (6.  5)  with  rQ  =  0  .  The  final 
phase  of  defoliation  occurs  after  both  traveling  hinge  circles  have 
reached  their  termTnafpo^itions  and  the  mechanism,  called  mechanism 
1,  is  described  by  (6.  5)  with  both  rQ  =  0  and  =  R  .  This  last  mode 
is  the  same  as  that  for  static  collapse. 

{b)  Mechanism  2a 

This  mechanism  given  by  the  velocity  field  (6.  5)  with 
r_  =  R  is  an  initial  mechanism  for  certain  ranges  of  a  and  X  .  In 
fact,  the  range  of  a  turns  out  to  be  0.  56  *  a  *  1  .  The  range  of  X 
depends  on  a  and  is  denoted  by  Xj  <  X  s  X2  the  numerical  values 
for  which  are  shown  in  Fig.  6.  5. 

(c)  Mechanism  2b 

This  mechanism  given  by  (6.  5)  with  rQ  =  0  is  likewise  an 
initial  mechanism  for  certain  ranges  of  a  and  X  .  The  range  of  a 
is  0  <  a  £  0.56  while  the  range  of  X  depends  on  a  and  is  again 
denoted  by  Xj  <  X  s  X^  the  numerical  values  for  which  are  shown 
in  Fig.  6. 5 

(d)  Mechanism  1 

This  mechanism, given  by  (6.  5)  with  r  =0  and  r7  =  R, 

o  - 

is  an  initial  mechanism,  and  hence  for  blast  pulses  the  only  mechanism, 
for  all  a  and  for  a  certain  range  of  X  depending  on  a  .  The  range 
of  X  denoted  by  1  <  X  s  Xj  is  obtainable  from  Fig.  6.5. 

The  way  in  which  X^  and  X ^  are  found  is  described 
once  the  governing  equations  are  derived. 
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FIG.  6.5  DEFORMATION  MECHANISM  DIAGRAM 


Governing  Equations 


Independent  of  mechanisms  the  equation  of  motion  is 


where  thu  subscripts  r  and  t  denote  partial  differentiation  and  m 
is  the  mass  per  unit  area  of  plate.  Differentiating  the  velocity 
(6.  5)  provides  the  acceleration  to  be  substituted  in  (6.  6)  and  the  yield 
condition  (Figs.  6.  3  and  6.4)  allows  the  circumferential  component  of 


bending  moment  N  to  be  eliminated.  Performing  the  integration  in 
(6.6)  and  using  the  properties  M=0  in  0  s  r  s  rQ  ,  Ki{r  ^  =  0  , 

Mf1^)  -  -  Mq  and  =  0  lead'  to  the  following  four  equations 

governing  mechanism  3. 


V'  =  oR  /12M 


(6.7) 


v'n(VV  Tl)[n(6-8T1+3r,2)  +  2(^-?2)(3-3n+T12)]-  V5/1T12[r1(l-T1)(4-3r1) 
+  (5j-?2)(6-8ri  +3ri2)]  -  V52  t]3(4-3ti  )  -  Vn  r|2[  r(4  -  3t) )  +  2(?j  -  §2)(3-2n)] 


-2?ll  25,  , 

1(?,-?2  +  n)2 

a  <  p 

2?1  2 
«  §2  +  r,)2 

pl  < 

(6.8) 

v'fSj-^+nJUe  (?1  ?2)-3-2(51-§2)(3-3t1  +  3t12-t,3)]  -  V51/lr3e2(?1‘?2) 


-V52l  3e 


r 52)^3-t,2(i- n)(3-2nA -2(5r52)2(3- 
Je  <5i“52-(i-ri)U(51-§2)(3. 

*>  am  m  \  y 


6r|  +  6ti2 -  2t]3)] 


h  "(1-’1))  +  (51-§2)(3-6ti  +6r,2-2r13)  +  3(l-n)] 


„  2(?r?2)  ,  ,  , 

-Vn  [3e  -  3  -  2  (§J-  ?2)(3-3n2+2n3)-  6(?1-?2)2(1  -  n)2J 

a<»- 

|l(pRV/4M0)^l-2;2+2tn(l/o)-e  ?2/a2J  -<1+ 5,- ?2>le2?I(5,-S2 +n) 


Pj  <  a 
(6.9) 
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r  2(5-?)  2(?  -?  ) 

V'iSj-Sj  +  iDUe  L  -3  +  2n(3-3n  +p)]-V?J[3e 

2  %  2  2 (5 1  " 

-3  -  2^ -  52)(3  -  6ri  + 6ri  -2r|)  +  2*1(1-11  )(3-2ri)]  V=2l3e 

r  ^  7  ,  2(?r§2) 

|?-(51-52)-l  +  2n  ?  +3-6t)  +  4n  -2ti]  -  Vn'[3e  -3-6(§1-?2) 

Z  2  i  /  r  2^  1  ”  Pj2^ 

(l-r!r  +  2T1<;(3-2T1)]  -  V?2[3e 

2  ? 

=  [(pR2a2/2Mo)  -  l]e  1(§1-?2  +  n)Z  0  <  p2  <  1 

(6.  10) 


The  new  variables  that  .wve  been  introduced  in  the  derivation  of 
(6.  tt),  (6.  9)  and  (6.  10)  are  defined  by 


j  =  tn(l/pj) 

S2  =  in(l/o2) 

11  =  1  -  PQ/P1 

o 

ii 

►1 

o 
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P1  =  *VR 

P2  =  r2/R 

Primes  over  the  dependent  variables  sr  ,  § j  .  52  *  and  *1  denote 
differentiation  with  respect  to  the  variable  t'  where 

r'  -  12M  t/mR^  (6.11) 

o 

The  equations  governing  mechanisms  2a,  2b  and  1  are  readily 
obtained  from  (6.  7),  (6.8),  (6  9)  and  (6.  10)  as  outlined  below. 

(a)  Mechanism  2a 

The  cuter  hinge  circle  is  stationary  at  radius  r  =  r2  =  R 
and  the  piopexties  leading  to  the  governing  equations  are  M  =  0  in 
the  region  0  s  r  &  rQ  ,  Mfrj)  =  0  ,  and  M(R)  =  -  Mq  .  Consequently, 
the  equations  are  (6.  7),  (6.  8)  and  (6.  9)  with  ?2  =  *2  =  0  . 
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(b)  Mechanism  2b 

There  is  no  inner  hinge  circle  and  the  relevant  properties 
are  M^)  =  0  ,  M(r2)  =  -  Mq  ,  and  Mf(r2)  =  0  .  Thus,  the  governing 
equations  are  (6.8),  (6.9)  and  (6.  10)  with  r)-  1  and  r|  -  0  . 


(c)  Mecha  ism  1 


The  outer  hinge  circle  is  stationary  at  radius  r  =  r2  =  R 
and  there  is  no  inner  hinge  circle,  leaving  the  two  properties  M(rj)  =  0 
and  M(R)  =  -  MQ  ,  so  the  governing  equations  are  (6.8)  and  (6.  9)  with 
?2  =  ?2  =  0  ,  ri=  i  ,  and  T)#  =  0  . 


Having  now  the  equations  governing  motion  in  all  four 
mechanisms,  it  is  possible  to  find  the  functions  Xj  and  \2  shown  in 
Fig.  6.  5.  It  is  assumed  that  the  hinge  circles  and  the  circle  r  -  r  j 

have  zero  initial  velocities  so  that  5(  =  5^  =  V  =  0  at  T  =  0  *  K1  is 
the  upper  bound  of  the  values  of  X  for  which  mechanism  1  applies  and 

the  initial  value  of  5j  »  designated  5j  ^  »  *or  e»ch  X  -  Pm/p§  or 
peak  pressure  pm  of  this  mechanism  is  obtained  by  solving  the 
following  two  equation  which  are  the  appropriate  special  cases  of 
(6.  8)  and  (6.  9): 


r 


V  '  (1+2^)  r.  / 


[  (P  R2 / 6M  )cx  (3-2ae  A)-l]e  (1+^) 
1  o  * 

7  *2?i  2?i 

[(p  R2/6Mje  -l]e  (l+?j) 
rm  °  1 


a  <  p, 


Pj  <  a 


(6.  12) 


2?1 

V/(3e  ■ 


f[(pmR2a2/2Mo)?1-(l+?1)le  Nl+Sj)  a  <P1 

_2?1/a2] 


•45  j)  =^[(pmR2a/4Mo)^l+2ui(l/a)- 


(6.13) 


2?, 


-  (1  +  ?!>)«  X<1  + 


Pj  <  a 
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By  considering  an  increasing  sequence  of  X's  the  upper 
bound  Xj  is  that  value  which  causes  an  inflection  point  in  the  bending 
moment  diagram  for  M  either  at  the  plate  center  or  at  the  support. 

This  means  that  the  yield  condition  will  be  violated  at  these  places 
for  values  of  X  infinitesimally  larger  than  X^  .  The  conditions 
Mrr(0,  0)  =  0  and  M^fR,  0)  =  0  are  expressed  by  the  equations 

v'  =  pmR2/12Mo  (6.14) 

,  2?l  1  ■>  2?i 

V  (3e  -1)  =  [pmR  a  /2Mo)- l]e  1(1+?J)  (6.15) 


and  when  X  =  X.  either  (6.  14)  or  (6.  15)  is  satisfied  by  the  solutions 
V '(0)  and  5^  of  (o.  12)  and  (6.  13).  Computations  show  that  when¬ 
ever  0  s  a  s  0.56,  X.  is  determined  by  the  condition  M  (0,0)  =  0 

i  r  r 

represented  by  (6.14),  and  whenever  0.  56  £  a  *  1  ,  X^  is  determined 
by  the  condition  Mfr(R,  0)  =  0  represented  by  (6.  15).  For  further 
increases  in  X^  the  former  case  indicates  that  the  initial  mechanism 
is  2b  whereas  for  the  latter  it  is  2a.  The  two  portions  of  the  Xj  curve 
are  shown  in  Fig.  6.5. 


The  values  of  X^  are  found  in  a  similar  way.  If  the  initial 

mechanism  is  2a,  the  initial  values  V'(0)  ,  ,  and  r|^  are 

obtained  by  solving  for  each  X(X>X^)  or  p  Equations  (6.  7),  (6.8) 

and  (6.  9)  with  the  special  values  =  ?2  =  11 '  ~  ®  1  an£*  ?2  =  °  '  ^2  *s 

that  value  of  X  which  causes  an  inflection  point  in  the  bending  moment 

diagram  for  M  at  the  support.  The  condition  M  (R,  0)  =  0  is 

r  r 

expressed  by  the  equation 


2? 


V'(3e  1-3  +  2ti(3 -3^ +n2)]  =  l(p 


R2C?/2M  )-l]e 
m  o 


2?, 


(1+Sj) 


(6.  16) 


and  when  \  =  initial  values  V  '(0)  ,  5^  and  satisfy 

(6.  16).  If  the  initial  mechanism  is  2b,  the  initial  values  V ' (0)  , 

(0)  1 

and  §2  are  obtained  by  solving  for  each  \  (\  >  \^)  or  pm  Equations 

(6.  3),  (6.  9)  and  (6.  10)  with  the  special  values  f  j  =  ^  =  h  '  =  0  ,  and 
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T1  =  i.  \  is  that  value  of  X  which  causes  an  inflection  point  in  the 

bending  moment  diagram  for  M  at  the  plate  center.  The  conditi  >n 

M  (0.  0)  =  0  is  expressed  by  the  equation 
rr 

V '  =  p  R2/12M  (6.  17) 

rm  o 

and  when  X  =  X2  ,  the  initial  value  V'(0)  satisfies  (6.  17).  The  two 
portions  of  the  X2  curve  are  shown  in  Fig.  6.  5. 

4.  Solution  for  Rectangular  Pulse 

In  this  section  is  an  outline  of  the  method  of  solution  applicable 
to  rectangular  pulses  with  pressures  high  enough  to  cause  initial  defor¬ 
mation  by  mechanism  3.  For  the  other  initial  mechanisms  the  method 
is  similar  and  the  details  simpler. 

A  solution  is  obtained  if  it  is  assumed  that  while  the  constant 

pressure  is  being  applied  the  hinge  circles  and  the  circle  of  radius 

r  =  rt  remain  stationary.  This  phase  of  the  motion,  called  phase  la, 

involves  the  solution  of  Equations  (6.  7),  (6.  8),  (6.  9)  and  (6.  10)  with  the 

pressure  p  a  constant  (X  >X_)  and  with  the  special  values  5i  =  ??  = 

^  10)  10)  (0)  ,  . 

t) '  =  0  .  The  solutions  are  denoted  by  >  ?2  ’  and  ^  .  V  is 

give,  explicitly  by  (6.7).  Corresponding  to  these  values  some  initial 

values  of  the  radii  pQ  ,  P1  ,  and  p2  are  listed  in  Table  6.  1  for  two 

values  of  a  .  Note  that  as  X  increases  they  tend  to  the  radius  a  of 

the  loading. 

Let  the  pulse  end  at  time  t  =  tQ  ,  or  when  t  =  tq  ,  and  let  the 
velocity  and  deflection  of  the  plate  center  at  this  time  be  VQ  and  6  fj  . 
Then  from  Equation  (6.  7)  alone,  successive  integrations  give 


=  I/m 

(6.  18) 

=  I^/Zmp 

(6.19) 

where  I  =  pt  is  the  impulse  applied  per  unit  area, 
o 
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INITIAL  VALUES  OF  p  .  P  ,  ,  AND 
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When  t  >  t  ,  no  pressure  acts  on  the  plate  and,  according  to 
(6.  7),  the  velocity  of  the  plate  center  or  rather  the  central  region 
0  *  r  *  rQ  ,s  constant  during  .his  phase  of  motion,  called  phase  lb. 

It  is  evident  from  (6.  8),  (6.  9)  and  {6.  10)  with  V '  =  p  =  0  that  §  j  ,  , 

and  r)  can  no  longer  be  treated  as  constants. 

In  Equations  (6.  8),  (6.  9)  and  (6.  10)  set  V  '  =  0  ,  V  =  V  and, 

o 

for  convenience,  introduce  the  nondimensional  time  variable  t  defined 
by 

t  =  mR2T'/12M  =  mR*' V  t/12M  (6.20) 

o  o  o 

The  resulting  equations,  in  indicial  notation,  may  be  written  in  the 
form 

■£.'  A. .  =  B.  i,  j  =  1,2,3  (6.21) 

in  which  the  derivatives  are  with  respect  to  t,  and,  for  notational 
purposes  only,  ?  3  =  r\  . 

The  chosen  numerical  scheme  for  the  solution  of  (6.21)  commences 
by  solving  for  the  derivatives  (Equations  (6.21)  are  linear  in  the 
derivatives)  to  give 

?/  =  •  ?2  ’  ?3>  i  =  1.2.3  (6.22) 

From  (6.21)  the  variable  t  can  be  eliminated  resulting  in  the  two 
equations 


d51/d?3  =  g(5r  ?2,  53)  and  d?2/d?3  =  h(51,52,?3)  (6.23) 

The  initial  values  determined  in  phase  la  are  2^  .  From  its  initial 
value  an  incremental  change  is  made  in  one  of  the  variables  here 
chosen  to  be  §3  so  that  its  new  value  is  §3^  =  5^  +  3  •  Then, 
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from  (6.23),  the  corresponding  new  values  of  and  are 

S',1*  *  «‘°>  *  g|0>  *5, 

-  5«°»  ♦  h«°» «, 

»her,  g'°>  .  ^  .  g<°>  .  g««»,  and  h'°>  .  h«<°>  .  5<°>  .  with 

(1)  1  C  °  (1)  (1)  1  ^ 
f  k  a  if«  Ir  am  P  ’  *  f  Vi  a  naui  i;n  I  itA  e  r»  ’  •  2LHd  H  m  k  a  _ _ J 


the  values  5-  the  new  values  g'  and  h'  9  can  be  determined  and 
hence  can  be  calculated  with  the  next  increment  “  5^. 

This  procedure  is  continued  until  either  § ^  =  rj  =  1  or  =  0  .  In  the 
former  case  the  deformation  continues  by  mechanism  2b  whereas  in 
the  latter  case  it  continues  by  mechanism  2a.  To  find  the  value  of  Tj 
of  t  at  the  end  of  phase  lb  all  the  time  increments  are  summed.  The 
first  of  such  increments  is 

m  =  -  T<°>  =  (?<■'  -  4°»)/;'<°» 

in  which  =  f^  by  (6.22). 

Finally  the  deflection  6^  of  the  plate  center  at  the  end  of  phase  lb 


is  given  by 


61  -  ‘o'  Vo«‘l-'o»  =  l2R2«Tl-’o»/12mM„  <6-24> 

2 

in  which  6  is  determined  by  formula  (6.  19)  and  t  =  12M  /R  p  . 

o  o  o  r 

As  mentioned  above  the  next  phase  of  motion,  called  phase  lc, 

can  take  one  of  two  forms  depending  on  which  hinge  circle  reaches  itn 

terminal  position  first.  Only  the  case  of  the  inner  hinge  circle  of 

radius  r  (t)  reaching  the  center  first  will  be  described  since  the  other 
o 

case  is  covered  by  the  description  of  phase  lb  in  the  problem  of 
Section  V  the  only  difference  being  in  the  initial  values  and  r)  . 

(In  Section  V  ?  plays  the  role  of  •)  Numerical  calculations  show 
that  rQ(t)  =  0  before  ^(t)  =  R  whenever  0  £  a  s  0.56.  The  equations 
governing  phase  lc  for  deformation  by  mechanism  2b  are  (6.8),  (6.9) 


and  (6.  10)  with  the  special  values  p  =  0  ,  =  0  and  r\  =  1  .  For 

brevity,  the  indicial  notation  is  used  to  express  these  equations  in  the 
form 

?/  C..  =  D.  i,  j  =  1, 2,  3  (6.25) 

i  iJ  J 


In  (6.25),  derivatives  are  with  respect  to  the  time  variable  t  defined 
by  (6.20)  and,  for  notational  purposes,  =  Y  =  V/Vq  . 


The  numerical  procedure  is  the  same  as  that  in  phase  lb  except 
that  increments  are  used  (instead  of  A5^)  stud  the  initial  values  are 

the  final  values  of  phase  lb.  The  procedure  is  halted  at  52  =  0  (p2  =  1) 
when  phase  lc  ends.  Let  the  central  deflection  and  the  time  at  the  end 
of  phase  lc  be  t>2  and  t2  (t  =  t2).  Then^ 

’2  -  T1  =X}T<n> 

n=  1 


where 


2T<n)  .  (?<"> 


-  5 


(n-l)WE'(n-l) 
2  '  ’2 


Also 


V6. 


Vdt  = 


iV  f 

i2mu;  j 


Y(t)ot 


which,  by  the  numerical  procedure  adopted  above,  is  replaced  by  the 
summation  jsj 


-  I MrT 

n=\ 


V<n,4t(n> 


(6.26) 


In  the  final  phase  of  motion,  called  phase  2,  deformation  is  by 
mechanism  1  with  the  governing  Equations  (6.  8)  and  (6.9)  having  the 
special  values  ?2=ri'=?2  =  p  =  °,  and  rj  =  1  .  So  that  the  analysis 
of  this  phase  conforms  to  those  of  Section  V  and  Reference  l  5.  l),  the 
variable  T  is  retained  and  the  variable  C  =  V/V2  is  introduced  where 
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=  Then,  ^e  iniii**  and  final  values  of  Q  are  unity  and  zero. 

Also,  the  governing  equations  become 


2  i 

C/(?1+D(251+l)  -  CSi?!  =  -(§!  +  !)  e  /v 


(6.27) 


2?i  2  3  ^5j 

C  (5j+l)(3e  -  3 -45j)  -  C§{  (3e  -3-6^  -2?p  *  -(?!  +  !)«  hz 

(6.28) 

in  which  differentiation  is  with  respect  to  t  and  y^  =  V2/Vo  *s 

the  value  of  y  at  the  end  of  phase  lc).  Eliminating  ^  from  (6.27)  and 

(6.28)  leads  to  a  linear  differential  equation  with  the  solution 


£.[  V^-J.xpl-  I  ...  ,J 

?l'T2)  +  1  J  4  +  75+25  -3e2^ 


(6.29) 


?l(T2) 


Motion  ceases  when  £  s  0  and  this  occurs  when  5]  «  0.478  which  is 
the  solution  of  the  equation 

4  +  75  +  252  -  3e2?=  0 


Let  t  =  Tj  when  motion  ceases.  Then 


VV 


t3  *  T  Z 


J  s'  2  J  «  *(?+ i)(4+7s+; 


-3)Cd5 


25“  -  3e^)  (6.30) 

5 1  (^2>  ^72* 

Finally,  let  the  central  deflection  when  t  =  t3  (t  =  t^)  be  6^  .  Then 


63_62  =  TOT  f  vdT  =  TItSt  Y2 
o  J  o 


5i(t3) 


f  (3e2?-45 
J  e2?(5+l)(4+ 


z-65-3)C2d5 


(5+l)(4+75+252-3ei^) 


W 


(6.31) 


In  Figs.  6.  6  and  6.  7  are  shown  the  paths  ABCD  followed  by  the 
point  (5j  ,  5^  ,  r()  for  values  of  a  equal  to  0.438  and  0.  656  respectively, 
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both  for  X  *  15  which  start*  motion  by  mechanism  3  (see  Fig.  6.  5  or 
Table  6.  1).  In  the  case  of  the  smaller  value  of  a  Fig.  6.6  shows  the 
trajectory  starting  at  A  and  intersecting  the  plane  rt  =  1  at  B  which 
corresponds  to  the  value  pQ  =  0  .  At  this  point  of  intersection  the  inner 
hinge  circle  has  reached  the  center  and  the  trajectory  continues  along 
BC  in  the  plane  q  =  1  which  corresponds  to  mechanism  2b.  It  next 
intersects  the  plane  §2  =  0  at  c  which  corresponds  to  the  value  p2  =  0. 
At  this  point  of  intersection  the  outer  hinge  circle  has  reached  the  support 
and  the  trajectory  continues  along  CD  ►  the  line  of  intersection  of  the 
planes  §-  *  0  and  T)  =  1  ,  corresponding  to  motion  by  mechanism  1. 
Figure  6.  7  for  a  =  0.  656  can  be  interpreted  in  a  similar  manner. 
Figures  6.8  and  6.9  show  the  trajectories  ABCD  in  Pq  ,  pj  ,  P2  space. 

Figures  6.  10  and  6.  11  show  the  variation  of  PQ  »  Pj  •  and  P2 
with  T  in  all  phases  for  values  of  a  equal  to  0.438  and  0.  656  respec¬ 
tively,  both  for  X  =  15  .  In  each  case  the  starting  values  of  Pj  are 
close  to  the  loading  radius  a  with  PQ  and  pj  almost  equally  spaced 
on  either  side  (see  Table  6.  1).  The  three  radii  are  constant  during 
phase  la.  The  magnitudes  of  the  average  velocities  in  phase  lb  are 
comparable.  During  phase  lb  and,  in  Fig.  6.  11,  during  phase  lc  the 
values  of  p'  are  almost  constant.  During  phase  2,  p  is  almost 
constant  in  Fig.  6.  1 1  for  a  =  0.  656  unlike  that  in  Fig.  6.  10  but  like 
that  in  Fig.  6.  12  which  has  been  included  to  allow  a  comparison  with 
the  case  a  =  1  (X  *  15.7).  The  final  value  of  t  ,  i.e.,  ,  increases 

with  a  which  is  to  be  expected  because  the  work  input  increases  with 
a  . 


Results  and  Conclusions 


For  three  values  of  a  Fig.  6.  13  chows  the  relationship  among 

the  final  central  deflection  6  ,  the  impulse  per  unit  area  I  ,  and  the 

2  2 

pressure  p  in  the  nondimensional  form  of  X  versus  mM  o6/rR“  = 
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FIG.  6.9  TRAJECTORY  OF  (p0,p,,p2)  FOR  o  -  0.656  AND  A  - 


FIG.  6.10  VARIATION  OF  pQlpv  AND  p2  WITH  r  FOR  a  -  0.438  AND  A  - 


FIG.  6.11  VARIATION  OF  p0,p ,,  AND  p2  WITH  r  FOR  o-  0.656  AND  A 
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FIG.  6.12  VARIATION  OF  p0lP],  AND  p2  WITH  r  FOR  a  -  1  AND  A  -  15.7 


The  curve  for  a  =  1.0  was  obtained  from  the  results  of  Section  V.  The 
three  curves  are  alike.  They  start  from  the  value  X.  =  1  where  6=0 
and  without  curvature  change  tend  monotonically  toward  vertical  asymp¬ 
totes.  The  location  of  these  asymptotes  have  not  been  found  for  a  *  0.438 
and  a  =  0.656  but  would  be  determined  by  considering  the  case  of  an 
ideal  impulse  (\  =  CO  ).  However,  judging  by  the  case  a  =  1.  0  for  which 
the  asymptote  is  known, the  value  of  v  when  \  =  100  is  sufficiently  close. 
These  values  are  0.  0212  and  0.  0456  for  a  equal  to  0.  438  and  0.  656 
respectively.  Figure  6.  14  is  another  way  of  representing  the  same  infor¬ 
mation  and  is  essentially  a  pressure -impulse  diagram. 

Figure  6.  15  is  a  nondimensional  plot  relating  the  pressure  to  the 
loaded  area  considering  as  a  parameter  v  which  is  proportional  to  the 
ratio  6/I2  .  By  treating  both  impulse  and  permanent  central  deflection 
as  fixed  quantities  the  curves  show  how  the  pressure  must  be  increased 
as  the  loaded  area  decreases  in  order  to  maintain  6  . 
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FIG.  6.13  RELATIONSHIP  AMONG  5,  A  AND  I 


From  Fig*.  6.  13,  6.  14  and  6.  15  the  following  conclusions  can 
be  drawn: 

(1)  For  a  given  impulse,  the  central  deflection  6  increases 
monotonically  with  the  pressure  Pxrx  becoming  a  maximum 
when  the  pressure  is  infinite,  that  is,  when  the  given 
impulse  is  applied  as  an  ideal  impulse 

(2)  Again  for  a  given  impulse,  consideration  of  how  6  increases 
with  \  in  Fig.  6.  13  reveals  that  for  a  >  0.438  the  central 
deflections  are  greater  than  35  percent  of  that  due  to  an  ideal 
impulse  for  pressures  about  sever  times  the  corresponding 
static  collapse  pressure,  that  is,  whenever  A  >  7  .  Below 
this  value  of  \  the  decrease  of  K  with  6  is  quite  pronounced. 
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(3)  For  a  given  central  deflection  6  ,  Fig.  6,  14  shows  that 
for  X  >  7  (a  >  0.438)  the  increase  in  impulse,  over  the 
ideal  impulse,  necessary  to  maintain  the  given  deflection 
is  less  than  7  percent,  A  large  increase  in  impulse  is 
required  to  maintain  6  as  X  decreases  further,  especially 
in  the  range  1  <  X  <  3  . 

(4)  For  a  given  permanent  central  deflection  6  and  a  given 
impulse  I  Fig.  6,  15  shows  that  the  pressure  required  to 
maintain  6  increases  rapidly  as  the  loaded  area  decreases, 
especially  for  a  <  0.  6  {for  the  values  of  v  shown). 


FIG.  6  *4  PRESSURE- IMPULSE  01  AC® AM 


6.  Description  of  Experiments 

Circular  plates  of  6C61-T651  aluminum  with  thicknesses  of  1/8, 
3/16  and  1/4  inch,  clamped  at  either  8-inches  or  12-inches  diameter, 
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were  subjected  to  pulses  generated 
in  an  oxyacetylene  shock  tube  having 
a  spring -mass  system  against  each 
plate.  The  experimental  arrange¬ 
ment  is  shown  in  Fig.  6.  16.  Thick 
steel  annular  supports  provided 
clamping  against  rotation  only,  a 
ring  spacer  keeping  their  distances 
apart  slightly  greater  than  the  plate 
thickness.  Apart  from  circumferen¬ 
tial  membrane  forces  and  unavoid¬ 
able  frictional  effects  the  plate 
material  was  allowed  to  displace 
freely  inwards  during  deformation. 
To  ensure  good  clamping  the  overall 
diameters  of  the  plates  tested  vere 
9  1/2  and  15  inches  but  to  minimize 
the  build-up  of  membrane  forces 
they  were  slotted  radially  at  about  1/2 -inch  intervals  around  the  rim. 

The  lengths  of  the  slots  were  such  that  their  ends  did  not  pass  over  the 
support  circle  during  deformation.  Figure  6,  17  shows  two  plates  of 
different  diameter  after  deformation. 

Two  types  of  pulse  were  used  and  their  shapes  were  taken  to  be 
those  recorded  under  a  similar  configuration  by  a  rod  gage.  They  are 
therefore  the  pressures  acting  on  a  fixed  target  and  they  imply  the  assump¬ 
tion  that  the  modification  of  the  pulse  due  to  plate  motion  is  small.  The 

2 

first  pulse  is  that  of  Fig.  2.8(d)  with  a  peak  pressure  p  =  865  lb/in 

m  2 

and  an  impulse  (area  under  pressure-time  curve)  1  *  0.  363  lb-sec/in6  . 
The  diagram  and  table  shown  opposite  the  pressure  record  gives  the 
configuration.  The  second  pulse  used  in  the  experiments  is  that  of 

Fig.  2.  9(e)  with  a  peak  pressure  p  *  485  lb/in  and  an  impulse 

•y  m 

I  =  0.209  lb-sec/in  . 


FIG.  6.15  VARIATION  OF  PRESSURE 
WITH  LOADED  AREA 
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FIG.  6.16  EXPERIMENTAL  ARRANGEMENT 


Tensile  tests  performed  with  the  Instron  machine  on  standard 
ASTM  specimens  taken  with  and  across  the  direction  of  rolling  of  the 
plate  material  provided  stress -strain  diagrams  from  which  the  yield 
stress  was  determined.  The  yield  stress  is  here  defined  as  the  stress 
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FIG.  6.17  TWO  PLATES  AFTER  DEFORMATION 


at  the  point  of  intersection  of  the  straight  line  approximations  to  the 
elastic  and  strain-hardening  portions  of  the  stress -strain  diagram. 

An  average  was  taken  of  the  yield  stresses  obtained  from  specimens 
with  and  across  the  rolling  direction. 

After  the  blast-loading  tests,  the  permanent  central  deflections 

6  were  measured, 
ex 

7.  Experimental  Results  and  Observations 

The  results  of  the  experiments  described  above  are  presented 

in  Table  6.  2  and  Fig.  6.  18.  Correlation  with  the  predictions  of  the 

rigid-plastic  theory  takes  place  through  the  ratio  &ex/&tk  °*  the 

central  deflections  in  Table  6.2  and  by  means  of  the  X  versus 

1/2 

IR/(12mM06)  plot  of  Fig.  6.  18  which  is  a  form  of  p-I  diagram. 

The  main  observation  is  that  the  central  deflection  ratios 
6 /6  ^  ,  by  ranging  from  0.041  to  0.416,  show  that  the  theory  over¬ 
estimates  the  central  deflections  from  the  present  experiments  by 
factors  ranging  from  about  2  1/2  to  24.  The  plates  are  much  stronger 
than  predicted  for  the  following  reasons: 

(1)  In  the  rigid -plastic  theory  the  effects  of  elastic  strain- 
energy  and  vibrations,  strain -hardening, and  strain-rate 
are  neglected,  all  of  which  add  strength  to  a  structure. 
Among  these  effects  it  is  probably  the  elastic  strain- 
energy  and  vibrations  that  are  dominant  in  the  present 
experiments.  The  plate  mate  rial,  606 1-T651  aluminum, 
was  chosen  because  it  shows  little  strain-hardening  and 
is  insensitive  to  strain-rate. 
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(2)  If  the  deflections  become  large  enough  membrane  forces 
become  significant  and  the  theory  considers  a  bending 
action  only.  However,  in  Table  6.  1  it  can  be  seen  that  the 
final  values  of  *.  / R  are  not  large  enough  to  allow  con¬ 
siderable  build -u^>  of  membrane  forces  l  1.  15J. 

(3)  Frict'onal  and  circumferential  membrane  forces  at  the 
support  give  the  plate  added  strength. 

(4)  The  experimental  pulse  is  not  rectangular  (see  Figs.  2.8(d) 
and  2.  9(e)).  For  a  given  impulse  a  rectangular  pulse  in¬ 
flicts  the  greatest  damage  among  blast  pulses.  This  has  not 
been  proved  in  general  but  is  true  in  all  known  solutions,  for 
example,  in  the  clamped  beam  problem  of  Section  3. 

(5)  The  applied  impulse  is  probably  less  than  that  recorded 
against  a  rigid  target.  Comparison  of  fixed  target  impulses 
obtained  from  rod^-gage  records  and  completely  free  or 
unrestrained  target  impulse  from  photodiode  velocity 
measurements  indicate  considerable  falling  off  of  impulse 
due  to  the  mobility  of  the  target  (see  Section  2).  It  is  also 
probable  that  the  pulse  shape  changes  with  the  mobility  of 

the  target. 


Among  the  above  reasons  for  the  unsatisfactory  correlation , the 
mobility  of  the  target  is  considered  the  most  important.  This  means 
that  the  impulse  imparted  is  less  than  that  assumed  and  that  the  deflection 
is  proportional  to  the  square  of  the  impulse.  For  ideal  impulses  on 
simply-supported  plates  of  6061-T651  aluminum  the  correlation  in  the 
range  of  6eJ£/R  used  here  is  &eJAh  -  6  [  1.  15]  whereas  the  best 

value  in  Table  6.  2  is  0.  42.  In  Table  6.  2  it  can  be  seen  that  for  a  given 
plate  radius  and  assumed  pulse  the  correlation  worsens  as  the  plate 
becomes  thinner.  This  observation  points  to  the  effect  of  target  mobility. 


The  interaction  of  the  elastic  vibration  modes  with  the  applied 
pressure  mechanisms  is  probably  also  very  important,  especially  since 
the  pulse  times  (see  Figs.  2.8(d)  and  2.  9(e))  are  greater  than  the  response 
time  or  quarter  period  T/4  of  the  fundamental  symmetric  mode  of  vibra¬ 
tion  for  each  plate.  For  the  plates  listed  in  Table  6.  2  the  response  times 

are  listed  in  Table  6.  3. 
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Plastic  moment 


NOMENCLATURE 


a 

d 

I 

m 

M.  N 


6 

ex 


loading  radius 
plate  thickness 
impulse 

mass  per  unit  area 

radial  and  circumferential  components 
of  moment 

fully  plastic  moment  0Qd^/4 

pressure 

peak  pressure 

static  collapse  pressure 

radial  coordinate 

radius  of  plastic  regime  boundary  (static 
solution) 

radii  of  plastic  regime  boundaries 
(dynamic  solution) 

plate  radius 

time 

pulse  duration 

time  at  end  of  phase  lb 

time  at  end  of  phase  lc 

time  at  end  of  phase  2  (motion  ceases) 

fundamental  period  of  plate  vibration 

velocity  of  plate  center 

V(tc) 

plate  deflection 
a/R 

w(o,  t~),  plate  central  deflection  at  time  t=t? 
(final  aeflection) 

w(o,  tQ) ,  plate  central  deflection  at  time  t  =  t 
w(o,  tj),  plate  central  deflection  at  time  t  =  tj 
w(o,  t^),  plate  central  deflection  at  time  t  = 
experimental  final  central  deflection 
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NOMENCLATURE  (Continued) 


T) 

X 


v 


Si 

h 

S3 


a 


o 


T 


T 


O’ 


/ 

T 


1  -  p0/p, 

mechanism  bounds  on  X 

mMo6/IZR2 

tn(l/ pA) 

£n(l/ p2) 

r\ 

r/P 

rQ/R  ,  Tj/R  ,  r2/R 

initial  values  of  pQ  ,  pj  ,  P2 
static  collapse  pressure 
[tn(r2/r1)  +  1  -  i^/rj]  -1 
yield  stress 
12Mot/mR2Vo 

12M  t  / mR2 V  ,  12M  t,/mR2V  .  12M  t,/mR2Vn  , 
00  o  ol  o  o  c  o 

12M  t./mR^V 
o  3  2  0 

12M  t/mR 

0  2 
12M  t  /mR 
o  o 
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SECTION  VII 


MATHEMATICAL  MODEL  OF  SPRING-MASS  SYSTEM 


1.  Introduction 

In  this  section  is  described  a  simple  mathematical  model  of  the 
layered  media  used  in  the  shock-tube  experiments  of  Section  2.  The 
experimental  configurations  and  corresponding  pressure -time  records 
are  shown  in  Figs.  2.  5  to  2.  9.  Among  these,  the  records  of  Figs.  2,  9(a), 
(d),  and  (3)  have  been  chosen  for  comparison  with  the  predictions  of  the 
mathematical  model.  In  the  photodiode  experiments  the  configuration 
of  Fig.  2.  9(e)  was  used  against  a  free  aluminum  disk  and  the  (x,  t)  plot 
of  the  disk,  shown  in  Fig.  .2.  16, is  compared  with  the  predictions  of  the  model. 

2.  Spring-Mass  System 

Instead  of  attempting  the  solution  of  the  complicated  problem  of 
finite  amplitude  wave  propagation  through  the  layered  media,  a  very 
simple  model  is  considered.  This  consists  of  a  spring-mass  system 
as  illustrated  in  Fig.  7.  1.  The  mass  of  foam  and  Mylar  are  assumed 

to  be  concentrated  at  the  ends  of 
each  massless  spring.  The  spring 
constants  are  based  on  the  stress- 
strain  relationship  of  polyurethane 
foam  (5  lb/ft^)  shown  in  Fig.  7.2. 
This  stress-strain  curve  is  itself 
approximate  since  it  was  obtained 
from  a  compression  test  on  an  un¬ 
confined  specimen. 

As  a  simplification  this 
stress-strain  curve  is  approxi¬ 
mated  by  two  straight  lines  or  a  straight  line  and  a  cubic  curve  as  shown 
in  Fig.  7.  3.  The  foam  is  assumed  to  compress  to  B  or  E  under  negligible 


F(t) 


Wj  in  ^ 


k,  k# 


U-4IU-H 


FIG.  7.1  SPRING-MASS-SYSTEM 
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FIG.  7.2  STRESS-STRAIN  DIAGRAM 

FOR  POLYURETHANE  FOAM 


FIG.  7.3  POLYURETHANE  SPRING 
PROPERTIES 


pressure  and  under  further  compression  behave  like  a  linear  (BC)  or 
cubic  (EF)  spring. 

If  the  displacements  of  the  masses  trij  ,  m^,  .  .  .  are  Xj  , 

.  ...  the  equations  of  motion  are 


in  which 


m  x  -  f  ,  +  f 
n  n  n-1  n 


n 


and 


n  n  ni 


k  (x  -x  , 
n'  n  n+r 


f  =  0 
o 


n  =  1 
n  0 

linear  spring 
cubic  spring 


(7.  1) 


The  dots  in  (7.  1)  denote  differentiation  with  respect  to  time.  If  there 
are  N  springs  and  N  masses  attached  to  a  rigid  target  (an  infinite 
N  +  1  th  mass)  the  final  equation  in  the  set  (7.  1)  consisting  of  N 
equations  requires 

/ 

linear  spring 
cubic  spring 


f 


n 
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If  there  are  N  springs  and  N  +  I  masses,  the  last  being  the 
relatively  hea  y  aluminum  disk  of  the  photodiode  experiments  the  final 
equation  in  the  set  (7.  1)  consisting  of  N  +  1  equations  requires 


Solution  of  the  Differential  Equation* 


In  the  set  of  Equations  (7.  I)  the  masses  m  are  obtained  from 

A 

the  densities  and  dimensions  of  the  styrofoam,  polyurethane  foam  and 
Mylar  (and  aluminum  disk  when  used).  The  spring  constant  k  or  £ 
depends  on  the  choice  made  for  B  or  E  in  Fig.  7.  4  when  fitting  the 
experimental  curve  in  Fig.  7.2.  It  also  depends  on  the  length  of  each 
foam  layer. 


The  force  applied  to  the  mass  is  taken  in  the  form 

F(t)  *  Fme"at  (7.2) 


The  numerical  value  of  Fm  and  a  are  based  on  the  pressure  caused 
by  a  shock  wave  from  an  oxyacetylene  gas  acting  on  a  rigid  wall  and 
were  obtained  by  fitting  curves  to  rod-gage  records  similar  to  that 
shown  in  Fig.  2.  3.  The  experimental  arrangement  for  these  measure¬ 
ments  is  that  of  Fig.  2.  1. 

The  initial  conditions  are  x  (o)  =  x  (o)  =  0  . 

n  n 

For  the  solution  of  Equations  (7. 1)  the  predictor  ard  corrected 
method  of  numerical  integration  has  been  used.  In  particular  the 
method  devised  by  Clippinger  [  7,  1 J  and  Hamming  [  7.  2]  has  been  used 
with  the  B5500  computer  [7.2,  7.4],  The  stability  of  the  methods  are 
discussed  in  References  [7.3]  and  [7.5]. 
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FIG.  7.4  EXPERIMENTAL  AND  MODEL  PULSES  —  CONFIGURATIONS  OF  FIG.  2.9(«) 


4.  Numerical  Results  and  Comparison*  with  Experiment 

The  first  comparison  to  be  made  is  with  the  pressure  record  of 
Fig.  2.9(e).  Corresponding  to  the  layered  media  of  the  experiment 
the  mathematical  model  had  10  masses  (and  an  eleventh  infinite  mass) 
and  10  springs.  Counting  from  the  loaded  end,  the  first  7  springs  were 
considered  cubic  with  £  =  26.  6  x  10^  lb/in'*  and  the  last  3  springs 
were  considered  linear  with  k  =  0.  014  x  10^  Ib/in,  the  latter  value 
being  taken  from  a  fit  of  the  lower  part  of  the  stress -strain  curve  of 
Fig.  7.2  because  the  pressures  are  lower  in  this  region.  An  attenuation 
factor  (reduction  of  the  peak  pressure  of  the  loading)  was  used  to  take 
account  of  the  plastic  nature  of  the  styrofoam  piston. 

Figure  7.4  shows  a  comparison  of  the  pressure  record  of 
Fig.  2.  9(e)  with  the  pressure  (on  the  fixed  plane  to  which  the  tenth 
spring  was  attached)  obtained  from  the  model.  Bearing  in  mind  the 
simplicity  of  the  model  the  agreement  is  quite  good,  at  least  for  times 
up  to  250  psec  from  the  arrival  time  of  the  pressure  wave.  After  this 
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time  agreement  deterioratei •  Certain  parameters,  such  as  the  value  of 
k  or  k  and  the  chosen  drop  in  pressure  to  simulate  the  effect  of  the 
styrofoam,  called  above  an  attenuation  factor,  were  adjusted  to  aid 
agreement  of  experiment  and  theory. 

A  similar  comparison  to  the  pulses  of  Figs.  2.9(a)  and  (d)  can 
be  seen  in  Figs.  7.  5  and  7.  6  and  again  agreement  is  quite  good  up  to 
250  psec  from  the  arrival  time  of  the  wave. 

In  the  photodiode  experiments  the  configuration  of  Fig.  2.9(e) 
was  used  against  an  unrestrained  disk.  By  means  of  the  photodiode 
technique  described  in  Section  II  the  (x,  t)  plot  of  the  disk  was  obtained 
and  is  showr  in  Fig.  2.  16  as  the  curve  labeled  "experiment".  The 
(x,  t)  plot  obtained  from  the  mathematical  model  is  also  shown  in 
Fig.  2.  16.  Since  the  transit  time  of  the  wave  was  shorter  in  the  case 
of  the  model  the  "experimental"  and  "mathematical  model"  curves 
were  positioned  by  giving  them  a  point  in  common  at  x  =  0.  0025  inch 
and  t  =  50  psec.  The  higher  transit  time  observed  in  the  experiments 
is  probably  mainly  due  to  the  friction  between  the  styrofoam  piston 
and  the  cylinder  wall  as  the  styrofoam  deforms  plastically.  Once  again 
it  can  be  seen  that  the  agreement  is  quite  good.  In  Fig.  7.  4  is  shown 
the  pressure  acting  on  the  disk  and  hence  the  predicted  effect  of  the 
mobility  of  the  target.  The  two  curves  labeled  "model  (fixed  target)" 
and  "model  (free  target)"  clearly  exhibit  a  divergence,  that  is,  the 
pressure  difference  increases  with  time  during  the  period  between 
t  =  200  psec  and  t  =  400  psec. 

It  is  concluded  that  the  simple  mathematical  model  gives  a 
reasonable  account  of  the  behavior  of  the  layered  media  used  in  the 
experimental  technique  for  pulse  shaping.  The  various  adjustments 
of  parametric  values  used  to  achieve  agreement  would  largely  become 
unnecessary  if  a  less  simple  but  more  accurate  model  were  constructed 
to  account  for  such  effects  as  wave  propagation  and  dissipation. 
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SECTION  VIII 


DEPENDENCE  OF  DAMAGE  ON  PULSE  SHAPES 


1.  Introduction 

In  this  section  a  comparison  is  rr.«»de  of  the  permanent  deforma¬ 
tions  of  simple  rigid-plastic  structures  caused  by  pulsek  of  equal  peak 
pressure  and  impulse.  It  is  proved  that  for  a  class  of  simple  rigid- 
plastic  structures  the  rectangular  pulse  causes  the  greatest  damage 
(defined  as  maximum  permanent  deformation).  Some  illustrative 
examples  are  given. 

2.  Simplest  Rigid -Plastic  System 

The  dependence  on  the  pulse  shape  of  the  maximum  displacement 
in  the  following  system  will  be  found.  Consider  a  pulse  p(t)  acting  on 
a  mass  m  per  unit  area  having  a  constant  resisting  pressure  p# 

(Fig.  8.  1).  Whenever  p(t)  becomes  larger  than  p#  the  mass  is  set 
in  motion,  this  motion  being  governed  by  the  equation 

p(t)  -  pg  =  mx  (8.1) 

where  x  is  the  displacement  from  the  initial  at  rest  positian.  Dots 
denote  differentiation  with  respect  to  time. 


FIG.  8.1  SIMPLEST  RIGID-PLASTIC  SYSTEM 
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With  the  initial  conditions  x(0)  =  x(0)  =  0  successive  integrations 
of  (8.  1)  yield 


I(t)  -  Pgt  =  mx 

(8.2) 

A(t)  -  p  t2/2  =  mx 
s 

(8  3) 

In  (8.  2),  I(t)  is  the  impulse  at  time  t  (area  under  pressure-time 
diagram  at  time  t)  and  in  (8.  3),  A(t)  is  the  area  under  the  impulse¬ 
time  diagram  at  time  t  .  Although  it  is  not  necessary,  it  is  convenient 
to  consider  pulses  with  an  initial  value  greater  than  p  ,  i.e.  ,  p(0)  >  p 

Let  the  mass  come  to  rest  at  time  t  =  t^  .  Then,  (b.  2)  with 
i(tf)  *  0  gives  t{  =  If/Ps  where  If  =  I(tf).  Substituting  this  value  of 
tf  in  (8.  3)  gives  the  final  displacement  in  the  form 

mXj  =  Aj  -  pgt2/2  =  Aj  -  Iftj/2  (8.4) 


in  which  A£  =  A(t^)  . 


By  means  of  expression  (8.  4)  the  values  of  x^  from  pulses  of 

equal  peak  pressure  pm  and  impulse  IQ  are  compared  to  that  from 

a  rectangular  pulse  of  pressure  pm  and  impulse  IQ 

difference  between  I  and  If  ,  they  are  deimed  by 

°  1  t 

CO  lf 


To  clarify  the 


•/ 


p(t)dt  and 


p(t)dt 


Two  cases  immediately  arise  depending  on  whether  Iq  -  1^  or 


I  >  L 


(1)  Case  1 


I  =  l< 
o  f 


In  this  case  the  whole  of  the  pulse  is  applied  before  motion 
ceases.  If  the  pulse  ends  at  time  t  =  t  ,  then  t  <  t.  .  Since 
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_ _ _ 


The  equality  in  (8.  5)  holds  only  for  a  rectangular  pulse  and  therefore 
is  a  maximum  for  a  rectangular  pulse.  It  was  shown  earlier  that  the 
term  Ijt^/2  is  the  same  for  all  pulses  of  this  case.  Hence,  from  (8.4) 
the  maximum  permanent  displacement  occurs  when  the  pulse  is  rectan- 
gular. 
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This  result  may  be  illustrated 

in  the  impulse-time  plane  of  Fig.  8.2. 

Af  for  a  rectangular  pulse  is  the  area 

under  whereas  for  all  other 

pulses  Af  is  the  area  under  the 

curve  OOj  and  line  OjFj  .  The 

triangular  area  under  OFj  is  1^/2. 

Thus  the  final  displacement  is 

proportional  to  the  difference  of  these 

areas  and  is  shown  shaded  in  Fig.  8.2. 

The  maximum  slope  of  the  curve  OO^ 

is  that  of  OO',  i.e.  ,  p  ,  and  OO, 

m  1 

lies  wholly  in  the  triangle  OO'Fj  the  area  of  which  represents  x^  when 
the  pulse  is  rectangular.  Note  that  the  slope  of  OF  j  is  pg  .  If  the  curve 
OOj  intersects  OFj  the  mass  comes  to  rest  (see  case  2  below). 

(2)  Case  2  I{  <  I 

In  this  case  pressure  is  still  being  applied  when  motion 

ceases.  Since  tr  =  I,/p  <  I  /p  ,  the  values  of  t ,  ,  unlike  case  1, 
f  i  s  o  s  f 

depend  on  the  pulse  and  are  all  less  than  the  value  of  t^  in  case  1. 
However,  the  reasoning  of  case  1  applies  leading  to  the  rectangular 
pulse  giving  the  maximum  displacement.  Expression  (8.  5)  becomes 
an  inequality  for  1^  <  I  and  A^  -  pgt^  /2  in  (8.3)  is  less  than  the 
value  corresponding  to  a  rectangular  pulse  at  the  same  time  t  =  t^  . 


FIG.  8.2  IMPULSE-TIME  DIAGRAM  — 
CASE  1 


An  illustration  of  this  result  can  be  seen  in  the  impulse - 
time  diagram  of  Fig.  8.3.  Since  I^/t^  =  pg  the  point  F£  lies  on  the 
line  OFj  of  Fig.  8.2.  The  area  under  the  curve  OF ^  is  A^  and  the 
triangular  area  under  the  straight  line  OF£  is  Ijt^/2  .  Their  difference, 
shown  shaded  in  Fig.  8.  3,  is  proportional  to  the  displacement  x^  .  The 
shaded  area  is  less  than  the  area  of  triangle  OO'Fj  which  is  proportional 
to  the  displacement  x^  from  a  rectangular  pulse. 


_ 


_ 
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Froi.i  the  above,  the  following 
theorem  may  be  stated.  Theorem: 
Among  all  pulses  of  equal  peak  pres¬ 
sures  and  impulse  th>e  rectangular 
pulse  causes  the  maximum  perma¬ 
nent  deformation  of  a  rigid-plastic 
structure  that  is  representable  by  a 
mass  with  a  constant  resisting  force. 

The  above  proof  can  readily  be 
modified  to  include  certain  pathologi¬ 
cal  pulses  such  as  double  pulses.  It 
is  believed  that  the  theorem  can  be  generalized  to  include  rigid-plastic 
structures  with  more  complicated  representations  or  modes  of  defor¬ 
mation.  This  is  illustrated  later  fox  the  case  of  a  beam  subjected  to 
uniform  blast  loading. 

A  few  examples  illustrating  the  applicability  of  the  above 
theorem  follow  immediately. 

3.  Rigid-Plastic  Cylinder 

Assuming  that  the  cylinder  remains  stable  [  1.9]  the  equation  of 
motion  when  a  pulse  p(t)  is  applied  uniformly  around  the  circumference 
is 


FIG.  8.3  IMPULSE-TIME  DIAGRAM  — 
TASE  2 


p(t)  -  0Qh/a  =  mw  (8.6) 

where  0q  is  the  yield  stress,  m  the  mass  per  unit  length  of  circum¬ 
ference,  h  the  thickness,  a  the  radius,  and  w  is  the  inward  dis¬ 
placement.  Dots  denote  differentiation  with  respect  to  time.  The  static 
collapse  pressure  is  pg  =  0Qh/a  and  is  regarded  as  the  constant 
resisting  pressure,  so  (8.  7)  may  be  written  in  the  form 

p(t)  -  pg  =  mw  (8.  7) 
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Two  integration!  of  (8.  7)  then  give 

mWf  S  Af  "  Pg4;  /2  (8-8) 

and  (8.8)  is  the  same  result  as  (8.  3). 

Consequently,  the  maximum  permanent  deformation  occurs  when 
the  pulse  is  rectangular. 

4.  Simply  Supported  Circular  Plate 

A  simply  supported  circular  plate  of  radius  R  is  subjected  to 
a  pulse  p(t)  uniformly  distributed  over  a  central  circular  area  of 
radius  a  .  This  is  the  problem  of  Section  IV  and  from  Eq.  (4.  10) 
the  final  central  deflection  is  given  by 

mw£  =  (12M0/pgR2)(Af-pgt£/2)  (8.9) 

Equation  (8.  9)  is  the  same  as  (8.  3)  so  that  the  rectangular  pulse 
causes  the  greatest  central  deflection  t p  is  given  by  (4.  6)). 

5.  Rigid-Plastic  Beams 

Consider  the  clamped  beam  problem  of  Section  III  in  which  a  blast 
pulse  is  applied  uniformly  over  a  central  part  of  the  span.  The  equation 
of  motion  governing  deformation  by  mechanism  1  is 

myQ  =3?o(l-?o/2)[p(t)-p8]  (8.10) 

in  which  yQ  is  the  central  deflection,  ?Q  =  a/f,  ,  where  a  is  the 
loaded  length  of  the  half -span  K  (see  Fig.  3.  1),  m  is  the  mass  per 
unit  length,  and  p#  is  the  static  collapse  pressure  given  by  formula 
(3.4) ,  and  the  dots  denote  differentiation  with  respect  to  time. 
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Two  integrations  of  (8.  10)  givs  the  final  central  deflection  in 
the  form 

my£  =  3?o(1-50/2)(A£  -  P(Jt*/2)  (8.11) 

Equations  (8.  10)  and  (8.  11)  are  similar  to  (8.  1)  and  (8.  3),  and  by 
applying  the  theorem  it  can  be  stated  that  a  rectangular  pulse  causes 
the  greatest  damage.  This  result  is  confirmed  by  the  results  of 
Section  III  in  which  the  central  deformations  due  to  a  rectangular,  trian¬ 
gular,  and  exponential  pulses  are  compared. 

6.  Rigid-Plastic  Beam  with  Moving  Hinges 

As  an  example  of  a  structure  which  does  not  quite  have  the 
representation  or  action  called  for  in  the  theorem,  a  simply  supported 
beam  subjected  to  a  blast  pulse  uniformly  distributed  along  the  entire 
span  is  considered.  It  is  shown  that  among  all  pulses  with  the  same 
peak  pressure  and  impulse  the  rectangular  pulse  causes  the  greatest 
damage.  This  is  an  indication  that  the  theorem  can  be  generalized. 

The  deformations  for  this  problem  have  been  found  by  Symonds 
[8.  l)  and  when  not  derived  the  required  results  will  be  extracted  from 
this  reference. 

For  peak  pressures  (assumed  to  occur  immediately  at  time 
t  =  0  (Fig.  5.  1))  greater  than  some  value  p^  to  be  determined,  the 
beam  deforms  by  mechanism  2  in  which  two  hinges  form  immediately 
and  travel  towards  midspan.  From  the  equation  of  motion  and  the 
continuity  of  velocity  condition  it  can  be  shown  that  the  location  of  the 
moving  hinges  is  given  by 


zo  =  6M0t/r  (8.  12> 

where  zq  =  L-xq  (Fig.  8.4),  Mq  is  the  fully  plastic  moment,  and  I 
is  the  impulse  per  unit  length  at  time  t  . 
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FIG.  8.4  BEAM  PROBLEM 


From  (8.  12)  the  initial  position  of  the 
hinge  in  0  st  x  <  L  is  given  by 

*2(o)  =  6M  /p  (8.13) 

o'  ’  o  m  ' 

From  (8.  13)  the  smallest  value  of 
Pm  for  which  traveling  hinges  occur  is 
obtained  by  setting  zq(o)  =  L  .  Let  this 
lower  bound  of  pm  be  pj  .  Then 

Pl  =  6M0/L2  =  3pg  (8.  14) 


where  pg  is  the  static  collapse  pressure. 

Whenever  p  *  p.  the  beam  deforms  in  mechanism  1  with  a 
m  i 

stationary  hinge  at  the  center  and  the  theorem  can  be  applied  directly. 

Consideration  is  given  here  to  the  cases  p  >  p,  because  the  theorem 

mi 

cannot  be  applied  directly. 

Let  the  moving  hingos  meet  at  midspan  at  time  t  =  tj  .  Then 
from  (8.  12)  and  (8.  14) 


tj  =  I/Pl  (8.15) 

Two  cases  are  now  considered.  In  the  first  the  pulse  ends  at 

a  time  t  =  t  which  is  less  than  the  time  t  =  t,,  that  is,  the  whole 
o  l 

of  the  pulse  has  been  delivered  in  phase  1 ,  while  the  beam  is  deforming 
by  mechanism  2.  In  the  second  the  pulse  extends  into  phase  2.  This 
second  case  may  be  subdivided  into  two  cases  depending  on  whether  or 
not  the  pulse  has  been  delivered  before  motion  ceases. 

(1)  Case  1:  Pulse  ending  in  phase  1,  tQ  <  tj 

In  this  case  the  central  deflection  and  central  velocity  at 
time  t  =  1 1  are 

Yo^i) s  Ai/m  and  y0(ti>  =  lolm 
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where  Iq  is  the  total  impulse  and  Aj  is  the  area  under  the  impulse - 
time  diagram  at  time  t  =  tj  .  Since  I  is  the  same  for  all  pulses  and 
Aj  is  a  maximum  for  rectangular  pulses  (see  proof  of  theorem)  phase  2 
has  with  the  same  initial  velocity  for  all  pulses  but  has  a  maximum 
Initial  deflection  when  the  pulse  is  rectangular.  Hence  the  rectangular 
pulse  gives  the  maximum  final  deflection. 

The  expression  giving  the  final  displacement  is  readily 

found  to  be 

myo(t2)  =  (3A2- Aj)/2  -  Io(3t2  -tj)/4  (8.  16) 

where  t2  =  3lQ/pj  =  3tj  is  the  time  when  motion  ceases. 

(2)  Case  2a:  Pulse  ending  in  phase  2,  t,  <  t_  <  t„ 

_  IOC 

During  the  time  tj  <  t  <  tQ  the  equation  of  motion  is 

mVQ  =  (3p-Pj)/2  =  3(p  -  ob)/2  (8.17) 

Noting  that  myQ(tj)  =  Ij  and  t^  =  Ij/p  from  (8.  15)  one  integration 
of  (8.  17)  gives 

mY0  =  (31  -  Pjt)/2  li  <  *  <  (®.  18) 

Noting  that  myQ(tj)  =  Aj  a  further  integration  gives 

myQ  =  (3A-Aj)/2  -  Pj(t2 -t])/4  tj  <  t  <  tQ  (8*  19) 

During  the  time  tQ  <  t  <  (t2  is  the  time  when  motion 

ceases)  the  equation  of  motion  is 

myQ  =  -  Pi/2  (8.20) 
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Making  use  of  the  continuity  of  velocity  and  displacement  at  time 
t  s  t  and  hence  making  use  of  (8.  18)  and  (8.  19)  successive  integrations 
of  (8.  20)  yield 

my0  =  (3Iq  -  Pjt)/2  tQ  <  t  <  t2  (8.21) 

myQ  *  (3A  -  Aj)/2  -  Pj(t2  -t2)/4  tQ  <  t  <  t2  <8-22> 

From  (8.21)  motion  ceases  at  time  t  =  given  by 

«2  =  3I0/P,  (8.23) 

and  from  (8.22)  and  (8.23)  the  displacement  at  time  t  =  is  given  by 

myo(t2)  =  <3A2  “  Al)/2  ‘  (3Io42  *  Vl*74  <8*24> 

Before  discussing  expressions  (8.  16)  and  (8.24)  that  for 
case  2b  will  be  derived. 

(3)  Case  2b:  Pulse  ending  after  motion  ceases,  t^  <  tQ 

During  the  whole  of  phase  2  the  equation  of  motion  is  (8.  17) 
and  the  initial  displacement  and  velocity  are  those  of  case  2a.  Hence 
the  expression  for  the  final  deflection  is.  from  (8.  19) 

myQ  =  (3A2  -  Aj)/2  -  pj(t2  -  t2)/4  (8.25) 


where,  from  (8.  18), 

42  =  3l27pl 


(8.26) 


Using  the  result  (8.  26)  for  t2  ,  (8.26)  may  be  written  in 

the  form 

myo  =  (3A2 -Aj)/2  -  (3I2t2  -  IjtjjM  (8.27) 

A  comparison  is  now  made  of  expressions  (8.  16),  (8.24) 
and  (8.27)  by  means  of  geometrical  considerations  in  the  impulse -time 
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planes  of  Fig.  8.  5(a),  (b)  and  (c).  For  this  purpose  it  is  convenient  to 
rearrange  (8.  16),  (8.24)  and  (8.27)  into  the  forms 


mVo(t2)  =  (A2  "Iot2/2)  +  1  (A2  "Iot2/Z)  •  {A1  -I0t1/2)]/2 
mVo(t2)  =  <A2  -W2*  +  I(A2’Iot2/2)  -  <A! -I^j/2)]  /2 
myo(t2)  =  <a2  ‘I2t2/2)  +  I  (a2  "I2t2/2)  '  (A1  -I1t1/2)]/2 


(8.24 


(8.27 


t0  *  m  p 


//  /  I 

/  /  1 
/  /  JX 
/  /  I 

I 


One  important  feature  of  the 
diagrams  is  that  the  lines  OR,  OP, 
and  OF  are  the  same  in  each.  OR 
is  the  path  that  is  taken  by  a  rectan¬ 
gular  pulse  and  has  a  slope  equal  to 
the  pressure  p  .  OP  has  a  slope 


(a)  CASE  I 


the  pressure  p  .  OP  has  a  slope 
2  m 

Pj  =  6Mq/L  and  OF  has  a  slope 
Pg  =  2Mq/L2  =  Pj/3  •  Three  typical 


r  " 

C0  —  7 - 

/  > 
"Tt( 

!//  i  . 


pulse  paths  OM  are  shown,  one  in 
each  of  cases  1,  2a,  and  2b. 


(b)  CASE  2 


By  algebraically  adding  the 
areas  represented  by  the  individual 
terms  in  each  of  (8.  16),  (8.24),  and 
(8.  27), it  can  be  seen  that  each  sum 
is  bounded  by  that  for  a  rectangular 
pulse,  this  sum  being  the  triangle 
ORF  plus  one -half  of  triangle  OPF. 
Thus  the  rectangular  pulse  causes 
the  greatest  deflection. 


(c)  CASE  3 


FIG.  8.5  IMPULSE-TIME  DIAGRAMS 
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radius  of  cylinder,  half  load  length  on  beam 
area  under  impulse -time  curve 
Aftj),  A(t2),  A(tf) 
thickness  of  cylinder 

impulse  (area  under  pressure-time  curve) 
total  impulse 
Ktj).  i(t2),  I(tf) 
half  span  of  beam 

mass 

fully  plastic  moment 

pressure 

peak  pressure 

static  collapse  pressure 

mechanism  bound  on  p _  for  beams 

radius  of  plate 
time 

pulse  duration 

time  at  end  of  phase  1 

time  at  end  of  phase  2 

time  when  motion  ceases 

duration  of  rectangular  pulse 

displacement  of  cylinder  or  plate 

displacement  at  time  t  =  t^  ,  w(t^) 

displacement  of  mass 

final  displacement  x(t^) 

coordinate  of  moving  hinge 

central  deflection  of  beam 

final  central  deflection  of  beam  y^t^) 

L-xo 

a/t  (a  is  half  load  length  on  beam) 
yield  stress 
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APPENDIX  A 


FORMULAS  FOR  THE  CENTRAL  DEFLECTION  OF  A 
CLAMPED  BEAM  DUE  TO  RECTANGULAR, 
TRIANGULAR  AND  EXPONENTIAL  PULSES 


In  the  following,  the  central  deflection  6  is  in  nondimensional 

form  and  denoted  by  v  where  v  =  6mM  /I?*2  •  The  notation  is  that 

o  l 

of  Section  HI. 


A.  Rectangular  Pulse 
(1 )  Mechanism  1 


v  =  3C2oa-50)2(K-l)/l6X 


(Z)  Mechanism  2a  and  3a 


v  =  ?0<2-50)(2K  -  1)/8X  -  1/12 


(3)  Mechanism  2b 

v  =  (§2/6)[(2-5o)/(2-?o-|oX)  +  2tn|(2X/3)/(2-50-^0XJ+21 

(4)  Mechanism  ?b 


>} 


v  =  (^/12)[5-4tn<250)-  3(2-So)/ZXSo] 


3.  Triangular  Pulse 
(1)  1  *  X  *  2 


0  s  C  *  1 
3o 


v  =  52(2-?o)2(X-1)3/X4 


<2)  2  s  X  «  X 


1 


0  s  5  s  1 
3  o 


v  =  52{2-§o)2(3X-4)/16X 
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(3)  s  U  2X  ?  2  1/2 

v  =  ?2,(2-?o)2(3X-4)/16\ 

f(2  -  ?  c)  (2  +  3?2  -  6?oK^  /  3  X )[  1-3?0(2-§o)/(6§o-3?2 _ 2)] 


(4)  X  2  2Xj  5q  *  1/2 

v  =  ?o(2-?o)(3X-2)/12X  -  1/2 


(5)  X  *  2X2  $o  £  1/2 

v  =  (?2/ 12)  [  5  -  4  tn  (2?q)]  -  ?q(2  -  ?q)/6X 


For  the  other  regions  of  the  failure  mechanism  diagram 


the  function  i  is  introduced. 

3 


where 


i  = 


c,V«l(b‘  -b‘)  +  64c  J  c2  tn(bj/b2) 
iil 


*1  *  32[6+  V/<2-?0» 

a2  =  -12  -  ?0l/<2-?0) 
a3  *  */3 

b i  *  8+  50X(2  -  t2)/(2-?0) 
b2  =  8  +  ;oM2  -T,)/(2-?o) 
c,  *  (8/3?2)[(2-?o)/x]4 
c2  =  8  +  25oW(2-50) 
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where 


(6)  s  \  s  2\j  0  £  5q  £  2/5 

Xj  *  X  *  \z  2/5  £  5q  £  1/2 

v  =  *  +  5q(2-50)2(3X  -4)/l6X  -  t2?2(2-?o)2(3X-3-TiX)/4X2 

Tj  =  2[1-3(2-?o)/X(2-3?o)] 
t2  =  0 

(7)  2Xj  s  X  s  X2  ?  £  2/5 

v  =  *  +  ?2/3  +  (52/3)tn[2X/(12  -650  +  35QX)] 

Tl=  1 
t2=  0 

(8)  X2  s  X  s  2X2  0  £  5o  £  2/5 

2Xj  £  X  £  2X2  2/5  £  5q  £  1/2 

v  =  *  +  (?2/3Rn[2X/(12  -  6?q  +  3?qX  )] 

+  ?ot2(2  ‘  5o)(3  -  t2)/6X  +  35o(2  ‘  ?0)2/16 

'  3Vo<2  -  5o)2(X  ‘  t1)/4x2 

where 

T,  -  M2  -  3?o)/6(2  -  5o) 
t2  =  2[  1  -  3(2  -  50>/S0\] 

(9)  \2  £  X  £  2X2  2/5  S  ?o  £  1/2 

v  =  $  +  §ot2(2  -  ?0)(3-t2)/6\  +?2(2-§o)2(3\-2)/8\ 

'  ?oTl<2  "50)Z(3  -  Tj)/4\-  3?2(2-5o)2(X2-4t2)/16X2 

where 

Tj  =  2[  1  -  3(2  -  ?C)A  (2  -  350)] 
r2  =  2[1  -3(2-50!/50X) 
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C.  Exponential  Pulse 


(1)  Mechanism  1 


v  =  (2  -  ?0)2(3;J/W:)[2M,£  +  e‘,£  -  1)  -  Tf2] 


where 


t,  =  X(l-e  *) 


(2)  Mechanism  2a  and  3a 


v  =  (2  -  ?0)(?o/4X)(Tl  +e  _1) 


+  (2-50)2{3?^/l6X2)[2X(Tf  +  e  'f-Tl-e  ' 1 )  -  Tf  +  if  ] 


where  T£  is  that  of  mechanism  1  and 

Tj  =  M65„  -3?2  -  2)(1  -«  T‘)/3;o(2  -  Co) 


(3)  Mechanism  2b 


v  =  (2  -  ?o)2  (?^/3) 


(l-e~T)2dT 


+  (2  -50)2  (352/16X2)[2X(Tf  +  e  'f-Tl-e  ' 1  )  -  jf  +  T2  ] 


where  t.  is  that  of  mechanism  1  and 
z 


Tj  =  X(2  -  350)(1  -e  1  )/3  (2  -  5) 


Mechanism  3b 


v  =  (2  -  sj«;/3) 


(1  -e~T  )2  dT 


+  5q  (2-?  q)(t2  +e  6  -  1)/4X  +  (2-?o)2(35  2/  16X2)[  2X(t£+  e 


■T1  "e 


1  .  2.2 


>-Tf +T1] 


l  I 


where 


1 


is  that  of  mechanism  1,  t,  is  that  of  mechanism  2b  and 
t2  ■  ?0xci  -  e  Tz)/3(2  -  ?o) 
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APPENDIX  B 


FORMULAS  FOR  THE  FINAL  DEFORMED  SHAPE  OF  A 
CLAMPED  BEAM  DUE  TO  A  RECTANGULAR  PULSE 


In  the  following,  the  deflection  y  is  in  nondime  ns  ional  form 
and  denoted  by  w  where  w  =  ymMQ/I^t^  and  §=x/(,  .  The  notation 
is  that  of  Section  III. 

A.  Mechanism  1 
w  =  3S^(1  -  ?)(2  -  ?o)2(X  -  1)/  16A 

B.  Mechanism  2a 

w  =  (1-5)  (2 +  5)/ 12-  (1-5g)2/4  -  S0(2-?o)/8\ 

™  =  (i-5)[(50/8X)(2-5o)/(i-52)  +  52/6  +  1/12] 

where  is  given  by 

U  -  ?2)2  -  3?c(2  -  ?0>/\+  3(1  -  ;o)2 

C.  Mechanism  Zb 

W  =  (5o/351)(51-5)[3(2-?o)/4\?1  +  1]  -(5*/3Rn§  0  s  §  *  5j 

w  =  (5^/3Rn(l/?)  5j  *5  *  1 

where 

§!  =  3 ?q/2  +  3(2  -  ?o)/2\ 


0 

h  *  5  *  1 
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D.  Mechanism  3a 


w 

=  (1  *  ?)(2  +  ?)/ 12  -(1-?0)2/4-?o(2-5o)/8X 

0  s  5  s  h 

w 

*  (l-?)(4?o-l)/12  +  (?0-?){2  -  350+  ?)/24 

+  <1-?o)2/24  +  ?0<2-50)/8X 

h  *5 s  h 

w 

=  (?0/8X)(2-?o)(?2-5)/(51-52) 

+  (1  -  ?)<4?0  -  ?j)/12  +  (1  -  ?Q)(1  -  §  j )/8 

h  *  5  ‘  Si 

w 

*  (l-?)(8?0-?-l)/24 

V) 

bS> 

V) 

where 

?,  =  «0  +  [3?0(2 -5^/X]* 

?2  =  «„  -  |3«0(*-S0»/Xl* 

52  *  2«„  -  * 

E.  Mechanism  3b 

*  =  (5^/12)[5-3)(2-5o)/2?oX-  tn(25o)] 

-  ?<2?q-?)/12  0  s  ?  s  ?2 

w  =  l50/8X)<2-5o)(?2-5,/(5l_52)+{?o/24)l9-3tn{2?o> 

+  2?1/?01  -  5(?2  +  250)/12  ?2*?* 

w  =  (50-?)(750-?)/24+  (?2/24)[5-8tn(25o)]  $j  *  ?  *  5j 

"  =  <?2/3Rn(l/5)  ^  *  ?  *  1 

where  5j  and  ?2  are  those  of  mechanism  3a  and  §j  =  2§q  . 
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